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SEISMIC  INTERFACE  WAVES  IN  COASTAL  WATERS:  A  REVIEW 

by 

Dieter  Rauch 


ABSTRACT 


The  basic  elements  of  seismic  wave  propagation  in  the  upper  sea  floor  are 
presented,  special  emphasis  being  devoted  to  the  often-neglected  interface 
phenomena  (Scholte-  and  Stoneley  waves).  The  most  important  theoretical 
approaches  to  the  seismic  sensing  of  waterborne  sound  are  discussed 
briefly  and  a  great  deal  of  the  pertinent  field  work  since  World  War  II  is 
reviewed.  Some  guidelines  are  given  for  future  research  activities  in  this 
interdisciplinary  field,  which  combines  the  low-frequency  regime  of 
underwater  acoustics  and  the  high-frequency  regime  of  seismology.  An 
extensive  bibliography  is  provided. 
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INTRODUCTION 


The  purpose  of  this  report  is  to  give  an  introduction  to  the  physics  of 
seismic  wave  propagation  in  the  sea-floor  and  to  present  the  state  of 
knowledge  on  seismic  sensing  of  waterborne  sound  according  to  the  open 
literature  published  before  1979.  Due  to  this  twofold  task,  the  report  has 
both  to  deal  with  different,  highly-specialized  fields  and  to  attempt  to 
synthesize  them  in  a  way  that  may  serve  as  a  starting  point  for  future 
interdisciplinary  studies.  As  a  result,  many  pertinent  aspects  of 
underwater  acoustics  and  seismology  can  be  discussed  only  in  a  quite 
superficial  way;  however,  the  interested  reader  is  furnished  with  an 
extensive  reference  list  to  compensate  for  these  unavoidable  restrictions. 

For  the  same  reason,  the  presentation  of  mathematical  derivations  is 
restricted  to  some  basic  but  quite  instructive  considerations  in  Sects.  1.2 
and  3.1.  The  reader  who  is  more  interested  in  the  phenomenological  or 
experimental  aspects  of  seismic  sensing  of  the  sea  floor  may  simply  skip 
those  chapters  without  interrupting  the  context. 

The  sound  field  radiated  by  a  water-borne  source  and  propagating  in  shallow 
water  or  from  deep  to  shallow  water  over  a  sloping  bottom  [i]  is  strongly 
affected  by  the  acoustic  properties  of  the  sea  floor.  The  ocean  bed  has  to 
be  considered  not  only  as  a  lossy  boundary  but  also  as  a  quite  complicated 
liquid/solid  structure  that  can  support  a  great  variety  of  wave-propagation 
phenomena.  From  this  seismic  point  of  view  it  is  obvious  that  the  acoustic 
energy  penetrating  into  the  sea  floor  may  sometimes  contribute  considerably 
to  medium-range  and  long-range  acoustic  transmission.  One  of  the  most 
drastic  examples  of  this  is  the  existence  of  a  low-frequency  cut-off  for 
the  first  mode  in  shallow  water.  For  all  frequencies  below  this  limit  the 
sound  field  ducted  by  the  water  column  is  in-phase  and  exponentially 
evanescent,  so  that  the  liquid  layer  loses  its  wave-guide  character.  The 
affect  of  a  very  unfavourable  sound-speed  profile  is  much  less  spectacular, 
of  course,  but  it  may  be  quite  obstructive  to  the  propagation  of  higher 
frequencies  as  well  and  thus  shade  off  shallow-water  areas  from  a  source  in 
deep  water. 

Under  these  restrictions,  even  medium-range  propagation  of  infrasonic 
etergy  (frequency  band  of  about  1  to  30  Hz)  can  happen  only  along  seismic 
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channels.  During  tne  enormous  progress  of  underwater  acoustics  since  World 
War  II  there  has  been  a  constant  but  rather  moderate  interest  in  these 
aspects  of  seismic  propagation  of  waterborne  sound.  Being  at  the  low- 
frequency  end  of  classical  sonar  activity  and  at  the  high-frequency  end  of 
seismic  research,  the  propagation  of  infrasonic  energy  was  tackled  only 
sporadically  and  incoherently.  Due  to  a  more  systematic  study  of  sea-floor 
acoustics  and  to  important  improvements  in  the  theoretical  means  (computer- 
models)  and  experimental  tools  (Ocean-Bottom  Seismometers  and  Seismic 
Arrays)  available  during  the  last  two  decades,  this  interdisciplinary  field 
seems  now  to  be  receiving  much  more  interest.  During  a  recent  workshop  at 
the  U.S.  Office  of  Naval  Research,  Washington  D.C.  [2]  the  experts 
complained  of  the  existing  lack  of  knowledge  and  reuummended  the  study  of 
two  potential  infrasonic  detection  methods:  either  the  deployment  of  ocean- 
bottom  seismometers  to  probe  the  marginal  portion  of  the  seismic  wave-field 
and/or  the  installation  of  rather  deeply  buried  seismic  arrays  (e.g.  in 
boreholes)  to  scan  its  shallow-refracted  portion. 

For  operational  reasons  (to  initiate  a  relatively  inexpensive  and  very 
flexible  system)  this  report  puts  emphasis  on  the  theoretical  and  practical 
aspects  of  the  former  project  (ocean-bottom  seismometers),  which  moreover 
fits  perfectly  into  SACLANTCEN ' s  current  sea-floor  studies.  Because  of 
their  efficient  excitatior  and  minimal  geometrical  spreading,  the  only 
seismic  wave-types  in  the.  marginal  zone  that  are  of  interest  are  those  that 
are  well  trapped  within  the  top  layers  of  the  ocean  bed  or  are  well-guided 
along  a  boundary  at  or  just  below  the  water/bottom  interface.  In  spite  of 
this  strong  guidance  the  waves  chosen  should  be  governed  by  the  least 
possible  dispersion  (minimized  multipath  propagation  or  almost  no  inter¬ 
action  with  other  nearby  boundaries)  and  no  radiation  losses  should  occur 
in  favour  of  other  wave  types.  This  report  will  show  that  apart  from 
horizontally-polarized  shear  waves  -  the  Love  waves  -  only  the  genuine 
interface  waves  of  the  Stoneley  or  the  Scholte  type  meet  all  these 
requirements  from  the  theoretical  point  of  view.  Accordingly  we  will 
discuss  their  physics,  examine  the  existing  experimental  data,  and  give 
some  modelling  results  by  using  SACLANTCEN's  Fast  Field  Program  (FFP). 

The  propagation-loss  curves  presented  are  calculated  with  realistic  (plane 
wave)  attenuation  coefficients  but  they  do  not  yet  account  for  parameter 
profiles  within  the  layers,  and  even  our  advanced  model  will  never  be  able 
to  take  into  account  spatial  variations  of  the  layers  or  inhomogeneities. 
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These  latter  features  alone  render  the  real-world  propagation  conditions 
much  more  complex  by  creating  a  favourable  upward  refraction,  and 
simultaneously,  very  detrimental  dispersion  and  scattering  effects.  Thus 
all  field  signals  will  strongly  depend  on  the  specific  terrain  traversed. 
Accordingly,  the  corresponding  receiver  configuration  (single  OBS-station, 
distributed  OBS-system,  or  OBS-arrays)  will  have  to  undergo  an  extensive 
calibration  procedure  before  it  can  be  used  for  surveillance  purposes. 
Concerning  the  sensors,  the  use  of  the  ocean  bottom  for  sound  transmission 
will  require  a  feasibility  study  of  seismometers,  alternately  called 
geophones,  instead  of,  or  together  with,  the  commonly-used  hydrophones. 

The  basic  characteristics  of  these  motion-detectors  will  therefore  be 
discussed  and  the  ambient-noise  background  in  the  sea  will  be  compared  with 
that  attributable  to  the  rather  unexplored  movements  of  the  sea-floor  in 
shallow-water  areas. 

Certainly  this  and  other  seismic  detection  methods  will  never  replace  the 
highly-developed  sonar  techniques  but  in  coastal  waters  they  may  offer  a 
valuable  support,  or  the  only  acoustic  alternative,  if  the  classical 
listening  techniques  lose  accuracy  or  even  fail. 
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1  THEORETICAL  ASPECTS  OF  SEISMIC  WAVE  PROPAGATION  IN  THE  SEA  FLOOR 

1 .1  Basic  Types  of  Seismic  Body  ar.d  Interface  Waves 

Due  to  its  latticed  microstructure,  a  solid  medium  can  support  two  basic 
types  of  elastic  waves:  the  congressional  or  P-wave  (sometimes  also  called 
"density"  wave)  having  a  particle  deflection  parallel  to  the  direction  of 
propagation,  and  the  shear  or  S-wave  having  a  particle  deflection  normal  to 
this  direction.  In  a  boundless,  non-absorbing,  homogeneous  and  isotropic 
solid  these  two  "body"  or  "bulk"  waves  propagate  independently  of  each 
other  with  their  characteristic  phase  velocities  of  Cp  and  c$  =  nCp,  with 
0  <  n  <  0.717.  These  fundamental  sound  speeds  (dealing  with  isotropic 
media,  we  will  understand  the  terms  "velocity"  and  "speed"  as  synonyms  and 
therefore  use  them  interchangeably)  are  related  to  the  mass-density  p  and 
to  the  elastic  constants  6  and  \i  (Lame's  constants)  or  a  and  E  (Poisson's 
and  Young's  modulus)  by 


(Eq.  la) 


(Eq.  lb) 


with  0  <  o  <  0.5  for  all  materials.  Another  often-used  constant  combining 


the  pair  (6,y)  or  (o,E)  is  the  modulus  of  compressibility 


K  =  3  .  3(l-2o)  (Eq.  2) 

SFPST  E 


Non-absorbing  liquids  have  no  shear  rigidity  (u=0)  and  can  accordingly 
support  no  shear  waves. 


Before  passing  over  to  more  complicated  wave  types  in  homogeneous, 
isotropic  media  we  would  like  to  clarify  some  often-used  notations.  In  the 
field  of  elastic  and  electromagnetic  wave  propagation  the  shortened  terms 
"plane  wave",  "cylindrical  wave",  or  "spherical  wave"  are  generally  used 
for  the  so-called  "homogeneous"  waves.  These  are  wave  fields  in  which  the 
surfaces  of  constant  phase  do  not  only  have  the  denoted  geometry  but  are 
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also  identical  to  the  surfaces  of  constant  amplitude.  Hence  we  will  in 
general  apply  the  same  nomenclature,  but  the  reader  should  keep  in  mind 
that  we  do  it  on  the  tacit  understanding  that  we  are  dealing  with  those 
homogeneous  waves  that  are  characterized  by  a  constant  amplitude  on  a  given 
phase  surface.  In  practice  almost  all  fields  are,  strictly  speaking,  "run- 
homogeneous",  having  surfaces  of  constant  phase  and  constant  amplitude  that 
do  not  coincide.  Accordingly  these  non-homogeneous  waves  are  provided  with 
a  more  or  less  irregular  distribution  of  the  amplitude  on  a  given  phase 
surface,  such  as  those  due  to  the  radiation  pattern  of  a  real  source,  for 
example.  It  has  to  be  remembered  that  the  present  report  is  essentially 
concerned  with  "inhomogeneous"  waves  representing  propagation  phenomena 
that  have  very  regular  and  well-defined  deviations  from  homogeneity. 

As  soon  as  "pure"  body  waves  impinge  on  a  free  sus  face  of  a  solid,  or  on  a 
marginal  interface  with  another  medium,  they  are  partly  converted  into  each 
other  and  their  coupling  may  create  "mixed"  wave  types  in  the  immediate 
vicinity  of  the  interface.  This  splitting  of  the  bulk  waves,  both  in  the 
reflected  ("inner  refraction")  and  in  the  refracted  field,  is  an  immediate 
consequence  of  the  specific  boundary  -  or  transition  -  conditions  that  have 
to  be  full f i 11 ed  there.  In  addition  to  these  conditions,  their  coupling  in 
the  form  of  different  interface  waves  requires  incoming  plane  or  cylindrical 
waves  with  complex  angles  of  incidence  or,  in  other  words,  the  grazing 
incidence  of  inhomogeneous  plane  or  cylindrical  waves  with  an  exponential 
amplitude  decay  for  growing  distance  from  the  interface.  At  first  glance, 
these  additional  requirements  seem  to  be  quite  extraordinary,  but  the 
mathematical  treatment  of  realistic  field  patterns  reveals  that  the 
occurrence  of  such  inhomogeneous  field  components  is  due  to  the  more 
complicated  (higher-degree)  curvature  of  the  incoming  homogeneous  wave- 
fronts  or  the  finite  width  (amplitude-shading)  of  the  incident  beams.  In 
the  simple  case  of  incoming  cylindrical  or  spherical  waves  we  can  see  this 
immediately  from  the  well-known  Sommerfeld  integral  [3  to  8]  or  Weyl 
integral  [5  to  9],  which  break  those  fields  down  into  an  infinite  sum  of 
inhomogeneous  plane  or  cylindrical  contributions. 

The  most  important  interface  phenomena  are  the  so-called  "free"  interface 
waves  that  propagate  along  the  surface  between  two  media  without  subsequent 
delivery  of  energy  from  sources  outside  the  interface.  We  have  already 
mentioned  that  their  existence  demands  the  combined  action  of  compressional 
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and  shear  waves,  which  is  why  at  least  one  of  the  two  media  has  to  be  a 
solid  while  the  other  may  be  a  solid,  a  liquid,  or  a  vacuum.  Accordingly, 
no  interface  wave  can  be  excited  at  the  surface  between  two  liquids,  except 
that  in  conditions  of  high  viscosity  there  may  occur  a  non-linear  propaga¬ 
tion  phenomenon  analogous  to  the  electro-magnetic  Zenneck  wave  [io]  at  the 
interface  of  two  dielectric  media. 

We  have  to  be  yet  more  precise  about  the  indispensable  shear  contribution 
to  the  synthesis  of  more  complicated  elastic  fields.  Only  that  component 
of  the  shear  deflection  that  is  polarized  vertically  to  the  given  boundary 
-  the  so-called  SV-wave  -  can  participate  in  the  formation  of  such  mixed 
wave  types  propagating  with  characteristic  trace  velocities  (=  resulting 
phase  velocities)  along  the  respective  wave  guide.  When  there  are  two 
adjacent,  homogeneous  half-spaces  (or  sufficient  high  frequencies), 
"genuine"  interface  waves  and/or  their  so-called  "generalized"  or  "pseudo" 
versions  may  exist.  The  genuine  interface  waves  have  a  constant  phase 
velocity  that  is  always  smaller  than  that  of  the  slowest  bulk  wave 
occurring  in  both  the  media,  and  their  amplitudes  decay  exponentially  on 
both  sides  of  the  common  surface  with  increasing  distance  from  it.  Thus 
they  represent  in  general  the  slowest  non-dispersive  field  component  that 
suffers  no  radiation  loss  in  favour  of  other  wave-types  and  restricts  its 
energy  flux  completely  to  the  immediate  vicinity  of  the  guiding  interface. 

Following  a  suggestion  of  Cagniard  [li]  the  genuine  interface  waves  go  by 
the  name  of  the  scientists  associated  with  them:  Rayleigh  [12] »  Scholte 
[13  to  15] ,  and  Stoneley  [16] .  These  are  listed  in  Table  1,  together  with 
an  estimate  of  their  respective  phase  velocities:. 


TASLE  1  CLASSIFICATION  OF  THE  BASIC  INTERFACE-WAVES,  WITH  AN 
ESTIMATION  OF  THEIR  PHASE-VELOCITIES 
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The  simplest  type  of  interface  wave  is  the  well-known  Rayleigh  wave,  which 
can  always  propagate  along  the  free  surface  of  a  solid  and  has  a  "penetration 
depth"  of  about  only  one  wavelength  [li,  17  to  21] .  Similarly,  there  are 
no  restrictions  on  the  existence  of  a  Scholte  wave  at  a  liquid/solid 
interface  and  its  decay  inside  the  solid  is  comparable  with  that  of  the 
Rayleigh  wave.  The  penetration  depth  inside  the  liquid  remains  equally 
small  if  the  adjacent  solid  is  acting  as  a  "very  slow"  or  "soft"  boundary 
(c'  «  c  ;  as,  for  example,  in  most  water/unconsolidated-sediment  combina- 

5  P 

tions),  but  it  can  be  much  bigger  if  the  solid  is  "extremely  fast"  or 
"hard"  (c*  >  c  ;  as,  for  example,  in  all  water/rock  combinations)  [6,  22  to 

s  p 

26].  The  most  complicated  type  of  interface  wave  is  the  well-known 
Stoneley  wave,  which  can  occur  at  the  common  surface  of  two  solids  for  only 
very  limited  combinations  of  parameters;  it  penetrates  into  each  solid  in  a 
similar  way  to  the  Rayleigh  wave  (27  to  31] . 

Concerning  the  mathematics,  all  these  genuine  interface  waves  result  from  a 
real  (or  almost  real)  pole  on  the  physically-relevant  sheet  of  the 
pertinent  Riemann  surface  ("top-sheet"),  which  means  that  their  velocities 
correspond  to  the  only  real  (or  almost  real)  root  of  the  pertinent 
equation.  For  a  liquid/solid  or  a  solid/solid  interface  there  may  be  some 
contribution  to  the  wave-field  from  an  additional  or  sometimes  vicarious 
complex  pole  on  a  physically  less  relevant  Riemann  sheet  ("lower  sheet") 
equivalent  to  an  additional  or  sometimes  unique  complex  root  of  the 
velocity  equation  with  physical  meaning.  This  special  singularity  can 
entail  a  generalized  or  pseudo-interface  wave  belonging  to  the  more 
extensive  class  of  "leaking"  interface  waves  [32] .  These  waves  are  so 
called  because  they  continuously  radiate  energy  into  the  "slower"  or 
"softer"  medium,  which  may  then  be  considered  as  a  kind  of  perturbation  of 
its  "faster"  or  "harder"  counterpart.  Because  of  their  radiation  loss 
these  wave  types  are  highly  attenuated  and  show  apparently  an  increasing 
amplitude/dis, ance  characteristic  on  the  side  of  the  "slower"  medium.  A 
well-known  example  is  the  pseudo-Rayleigh  wave  propagating  along  a 
liquid/solid  interface  with  a  phase  velocity  practically  identical  to  that 
of  the  genuine  Rayleigh  wave  [ 23 ;.  For  a  great  many  water/solid  combina¬ 
tions  the  existence  of  this  leaking  wave  has  been  demonstrated  with  the  aid 
of  seismic  records  or  tests  in  the  laboratory.  However,  if  the  liquid  is 
relatively  "fast"  and  the  adjacent  solid  not  "hard"  enough  (as,  for 
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example,  in  the  case  of  most  water/unconsol i dated-sediment  combinations) 
this  pseudo-Rayleigh  wave  is  completely  suppressed  in  favour  of  the  genuine 
Scholte  wave  (see  also  Sect.  1.5).  Much  less  known,  and  quite  unexplored, 
is  the  pseudo-Stoneley  wave  that  occurs  at  the  common  surface  of  two  solids 
and  might  be  an  alternative  to  overcome  the  restrictions  imposed  on  the 
existence  of  the  genuine  Stoneley  wave.  Unfortunately,  the  nomenclature 
concerning  these  interface  phenomena  is  still  not  yet  standardized  and  many 
authors  call  all  of  them  Rayleigh  waves  or  distinguish  merely  between 
Rayleigh  and  Stoneley  waves.  This  has  created  considerable  confusion  and 
misunderstandings  about  the  variety  of  wave  types  described  above. 

Finally,  we  should  discuss  briefly  the  neglected  component  of  the  shear 
deflection  that  is  polarized  parallel  to  a  given  interface  —  the  so  called 
SH-wave.  This  type  occupies  a  special  positiion  in  so  far  as  it  needs  no 
other  wave  to  fulfil  the  boundary  or  transition  conditions.  As  a  self- 
consistent  modification  of  the  shear  wave  it  creates  the  only  "pure" 
propagation  modes  (see  Sect.  1.4)  within  solid  layers,  which  are  named  Love 
waves.  Because  they  are  perfectly  reflected  or  guided  by  a  free  surface  or 
a  liquid/solid  interface  these  Love  waves  represent  an  important  factor  in 
seismic  wave  propagation  and  are  often  included  with  the  interface 
phenomena  [5.29,33,34] . 

To  give  examples  of  seismic  interface  waves  and  to  show  their  integration 
into  the  physical  background  we  will  refer  several  times  within  this  report 
to  the  medium  or  microscopic  scale  of  elastic,  electro-magnetic,  and 
optical  phenomena;  that  is,  to  the  very  high  frequency  regime.  Just 
recently  the  existence  of  piezoelectric  acoustic  boundary  waves  has  been 
proved  [35,36],  and  some  years  ago  even  more  elementary  analogues  have  been 
discovered  in  quantum  mechanics:  surface-spin  waves  in  semi-infinite 
ferromagnetic  crystals  [37]. 


1 .2  Elementary  Relations  for  an  idealized  Plane  Harmonic  Scholte  Wave* 

To  give  some  insight  into  the  physics  of  the  interface  problem  this  section 
presents  a  brief  mathematical  description  of  a  plane,  monochromatic  Scholte 


*Th<5  iee.ti.cm  ii  baied  on  a  SACLANTCEN  inteAnal  working  [xiper  UN-461)  that 
um  uriginatty  written  in  VecembeA  1975  and  tateA  r epAoduced  cu  SACIANTCEN 
Memorandum  SM-1 38  in  June  1980, 
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wave  propagating  along  the  adjacent  surface  of  two  homogenous,  isotropic, 
and  non-dissipative  half  spaces.  The  results  give  a  clear  idea  of  the 
pertinent  propagation  mechanism,  but  they  are  not  suitable  for  the 
treatment  of  a  realistic  sea-floor  environment  where  parameter  gradients, 
inhomogeneities,  stratifications,  and  material  losses,  for  instance, 
provoke  serious  modifications. 

We  choose  Cartesian  coordinates  (x,z)  for  range  and  depth  and  distinguish 
between  water  (z<0)  and  solid  bottom  (z>0)  by  the  use  of  the  subscript  "w" 
for  all  terms  relating  to  the  liquid  half-space.  If  we  consider  the 
simplest  disturbance  -  a  plane,  monochromatic  wave  of  angular  frequency 
10  =  2nf  propagating  in  the  +x  direction  -  our  problem  becomes  two- 
dimensional  (^-  =  0).  Therefore  the  vector  of  particle  displacement  is 
reduced  to  s  =  (^,0,w)  and  the  time-dependence  to  e-la)  .  (Here  and 
throughout  this  discussion,  all  underlined  symbols  represent  complex  terms 
according  to  the  definition:  a  =  a+i a '  and  derivatives  with  respect  to 
time  t  are,  as  usual,  denoted  by  a  dot). 

If  we  remain  in  the  framework  of  linear  elasticity,  neglect  superimposed 
body  forces,  and  introduce  the  displacement  potentials  4>  and  ip  according 
to: 

s  =  (£,0,w)  =  grad£  +  rotj  , 

with  (Eq.  3) 

|  =  (0,M), 


the  above-mentioned  potentials  have  to  fulfill 

water:  A<j>  +  k,2* ,  =  0  with  ktl  =  — 

-w  w—w  w  c, , 

w 

Aj>  +  k^  =0  with  kp  =  £- 

bottom:  P 

Ajj  +  k2^  =0  with  k  =  -p- 

s  s  cs 


the  wave  equations: 

(Eq.  4a) 

(Eq.  4b) 
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For  their  solutions  we  can  calculate  the  non-vanishing  components  of  the 
stress-tensor  X  according  to: 


water: 


t  =  T  5 

-xxw  "xxw 


T  =0 
-XZ 

w 


a, 

3x2  3z2 


(Eq.  5a) 


bottom: 


324>  a2<)) 


32<f>  a 


=  S( — -  +  — )  +  2y  ( - ) 

3x2  3z2  3x2  3x3z 

32(J>  32i|;  32i|< 

=  U(2  —  +  ---=) 

3x3z  3x2  3z2 

3  24>  32*  324>  d2ip 

=  6(— -  +  — -)  +  2p  (—  + - ) 


3x  3z 


3z2  3x3z 


(Eq.  5b) 


In  our  case  the  transition  conditions  for  the  sound  field  at  the  interface 
(three  boundary  conditions  for  three  force-free  potentials)  are  given  by: 


z  =  0 


w  =  w 
— w  — 


£  =  Izz 
Ixz=  0 


(Eq.  6) 


If  we  assume  displacement  potentials  of  the  form: 


water:  *  =  Me*2  ei(kx-“t> 

-4V 


(z  <  0) 


(Eq.  7a) 


*  =  ge'1*2  e’(kx-u>t) 
bottom:  *  =  Se"sz  ei<kx-“t> 


(z  »  0) 


(Eq.  7b) 


with  k  =  - 


to  fulfil  the  wave  equations  (Eq.  4),  the  vertical  wave  numbers  describing 
the  transverse  decay  of  the  field  have  to  be: 


k2-k2 ,  q  =  /k2-k2 ,  s  =  /  k2-k2 


(Eq.  8) 
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The  use  of  the  continuity  conditions  (Eq.  6)  leads  finally  to  a  set  of 


three  equations 

for  the  amplitude 

/  2ikq 

(k2+s2) 

0 

( u(k2+s2) 

-2i'hiks 

6w(k2-m 

^  q 

-ik 

m 

(Eq.  9) 


According  to  the  first  and  the  last  equation  these  amplitude  factors 
are  related  by: 


m  k2+s2 

q  1,2 

Ks 


M 


(Eq.  10) 


The  set  of  homogeneous,  linear  equations  (Eq.  9)  has  a  non-trivial  solution 
only  if  the  coefficient  determinant  is  vanishing  (in  the  case  of  more 
complex  systems  the  same  procedure  leads  to  the  pertinent  dispersion 
relations  from  which  the  characteristic  phase  and  group  velocities  are 
calculated.)  With  the  abbreviations: 


R  =  if) 

P 


N-  = 


l-2o 

2 


CS  ^ 

x  =  if) 


the  decay  constants  (Eq.  8)  take  the  form 


(Eq.  11) 


k.lyT 

m  =  k$  vO(-N‘,  q  =  kg  A-R,  s  =  kg  /X-V 


(Eq.  12) 


and  the  vanishing  determinant  leads  to  the  following  equation  for  X  and 
thus  to  the  phase  velocity  c  of  the  Scholte  wave 

/X^R- 

4X/X-1  »^R  -  (2X-1)2  =  H -  (Eq.  13) 


For  the  limiting  case  H  =  0  (i.e.  p, ,  =  0),  the  right-hand  side  disappears 

w 

and  we  end  up  with  the  well-known  Rayleigh  wave  equation,  which  has  always 
one  positive  real  root  c  =  cR  <  cs  (see,  for  example,  [6,18,22]). 
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The  latter  phase  velocity  is  frequency  independent  and  can  be  approximated 
to  high  accuracy  by  a  simple  formula  [38]: 


„  _  0.87  +  1.12a  _ 
CR  =  TTa  cs 


(Eq.  14) 


For  H  ^  0,  Eq.  13  describes  a  coupled  system  of  two  components,  each  of 
which  contains  "one  degree  of  freedom",  which  is  why  we  have  to  expect  two 
different  roots  with  physical  significance  [6,11,15,22] .  The  first  root  is 
always  complex  and  characteristic  for  the  so-called  pseudo-Rayleigh  wave 
with  the  phase  velocity  c  =  cR  +  icR;  the  second  is  always  positive  real 
and  belongs  to  the  Scholte  wave  c  =  cSc^.  (To  find  those  solutions 
numerically  it  makes  no  sense  to  square  Eq.  13  and  to  profit  from  standard 
subroutines  for  digital  computers,  because,  as  this  procedure  leads  to  a 
polynomial  of  the  eighth  degree  in  X,  it  may  be  extremely  difficult  to 
separate  the  two  pairs  of  "physical  roots"  for  the  genuine  and  the  pseudo¬ 
interface  from  the  other  six  pairs  of  extraneous  ones.) 


If  the  adjacent  fluid  has  a  very  low  density  (H«l,  such  as  air,  for 
example)  and  thus  represents  merely  a  "slight  disturbance"  of  tne  solid, 
only  the  pseudo-Rayleigh  wave  will  be  excited.  As  the  sound  speed  cQ  in 
such  a  gaseous  atmosphere  is  usually  smaller  than  the  phase  and  trace- 
velocity  cR  «  cR  of  this  interface  wave,  a  slight  amount  of  energy  is 
continuously  transferred  to  a  dragged  sound  field  (head  wave)  in  the  fluid 
medium.  Due  to  these  radiation  losses  the  pseudo-Rayleigh  wave  shows  an 
additional  (non-dissipative)  decay  of  the  amplitude  in  the  propagation 
direction.  This  decay  causes  a  fast  extinction  of  that  wave  type  in  most 
liquid/solid  systems  with  cQ  <  cR  (such  as  at  a  water/rock  interface), 
because  the  higher  density  of  typical  liquids  creates  a  much  stronger 
coupling  (better  impedance  matching)  and  thus  a  more  intense  energy  drain. 
As  a  conseqence  of  this  peculiarity,  the  pseudo-Rayleigh  wave  can  no  longer 
be  excited  if  the  sound  speed  of  the  liquid  reaches  or  exceeds  the  trace 
velocity,  i.e.  c0  ^  cR  (such  as  at  most  water/sediment  interfaces).  In 
this  connection  the  reader  may  think  of  the  famous  "air-coupled"  Rayleigh 
wave  often  generated  by  earthquakes  with  shallow  epicentres  or  by  big 
explosions  near  to  the  earth's  surface.  However,  this  most  prominent 
signal  of  many  seismic  records  is  a  coupling  phenomenon  of  two  independent 
wavefronts  and  occurs  merely  as  a  sort  of  "non-radiating"  limit  of  the 
pseudo-Rayleigh  wave.  This  coupling  of  the  shock  wave  in  air  with  a 
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"speed-synchronized"  wavelet  of  a  realistic  Rayleigh  wave  pulse  in  the 

bottom  is  due  to  the  unavoidable  dispersion  of  seismic  interface  waves  (see 

Sect.  1.4).  Thus  the  wavelet  of  the  "air-coupled"  Rayleigh  wave  separates 

the  pulse's  low-frequency  head,  with  pseudo  Rayleiyh-wave  character,  from 

its  high-frequency  tail,  with  almost  "pure"  Rayleigh-wave  features.  f 

When  the  liquid  has  a  lower  or  a  higher  sound  speed  than  the  Rayleigh 
velocity  of  the  solid,  the  Scholte  wave  becomes  the  more  important,  or  all- 
important,  propagation  mechanism.  As  a  more  detailed  investigation  has 
shown  in  [30]  the  resulting  phase  velocity  csc^  of  this  wave  type  is  also 
frequency-independent  (non-dispersive  propagation),  and  always  slightly 
smaller  than  the  lowest  velocity  occurring  in  both  media  i.e,  c$c^  <  Min(cQ,  cR) 

(see  Table  1).  The  interested  reader  may  find  many  normalized  solutions 
X4  or  ( X/N) ~ ^  of  Eq.  13  in  the  paper  of  Strick  and  Ginzbarg  f 24 ;,  where 
they  are  plotted  as  a  function  of  the  velocity  ratio  fT^,  with  the  Poisson- 
modulus  a  and  the  density  ratio  H  as  parameters.  Thus  velocity  and 
wavelength  of  the  Scholte  wave  amount  almost  to  the  corresponding  entities 
of  a  frequency-synchronized  water  wave,  if  we  deal  with  an  "extremely-hard" 
sea  floor  (c  >  c1(),  such  as  mother-rock,  for  example.  Normally,  with  a 
"relatively  hard"  or  a  "soft"  ocean  bottom  (c  <c.,  or  c,«c,,),  such  as 

i>  i%  D  W 

sedimentary  rocks  or  marine  sediments,  these  entities  are  much  smaller  and  4 

may  even  decrease  to  a  few  percent  of  those  of  the  water  wave.  With  the 

use  of  hydrophones  or  geophones  one  can  detect  the  omni-directional 

pressure  within  the  water  mass  jr  the  vector  of  particle  displacement  (or 

its  time  derivatives:  particle  velocity  and  particle  acceleration).  For 

this  reason  we  now  determine  the  real  components  of  the  displacement  vector 

and  of  the  stress  tensor  caused  by  a  plane,  monochromatic  Scholte  wave. 

With  the  aid  of  Eqs.  3,  7  and  12  we  find  the  following  displacement 
components: 

uw  =  %  sin(kx-a)t) 
with 

Gw  =  -ksMYemz 

water: 

A 

w,,  =  w,  cos(kx-wt) 
w  w  '  ' 

with 

Sw  =  ksN/ff  emz  =  ksQ  ^  e"z  >  0 


■k  0  ■  -=£I— 
S  Ah^(2X-1  ) 


,mz 


<  0 


(Eq.  15a) 
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u  =  u  sin(kx-u)t) 


with 


bottom: 

w  =  w  cos(kx-tot) 


u  =  k.QY  ^  e-«  .  e-iz  |  <  0 


(Eq.  15b) 


with 


w  =  k 


.d/wf  { ^  e~sz  -  e-£lz  }  >  0 


The  results  are  in  each  case  the  parametric  representation  of  an  ellipse 

having  its  main  axes  parallel  to  the  axes  of  our  coordinate  system.  With 

increasing  distance  from  the  interface,  the  displacement  amplitudes  u  .  w 

w  w 

and  w  decrease  exponentially  without  changing  sign.  The  horizontal 
amplitude  in  the  bottom,  u,  shows  the  same  asymptotic  behaviour  after 
having  changed  sign  at  a  depth  of  the  order  of  one  tenth  of  the  pertinent 
wavelength.  Figures  la  and  b  give  qualitative  sketches  of  these  particle 
displacements  and  their  resulting  orbits  for  a  typical  water/sediment 
interface.  While  the  "penetration  depth"  in  the  solid  is  always  limited  to 
a  range  or  about  one  wavelength  (similar  to  that  of  the  Rayleigh  wave),  it 
is  slightly  or  substantially  bigger  on  the  side  of  the  liquid  if  the  solid 
medium  is  acting  as  a  "relatively  hard"  or  a  "very  hard"  boundary 
(relatively-bad  or  extremely-bad  impedance  matching).  As  the  sea  floor  in 
coastal  waters  consists  mainly  of  unconsolidated  and  consolidated  sediments 
we  may  expect  that  waterborne  infrasound  is  transmitted  as  a  Scholte  wave 
even  in  very  shallow  water  areas. 


The  interesting  components  of  the  stress  tensor  can  be  derived  from  Eq.  5 
by  referring  once  more  to  Eqs.  7  and  12: 


P 

water: 


Pw  cos(kx-wt) 


with 


A 


■  -Vi*' 


mz 


pHv^R  gmz 
v^N  (2X-1) 


(Eq.  16a) 
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Txx  =  Txx  “sO«-"t> 

with  _ t 

Jxx  ‘  ksO  {  e"SZ  -t«IWM]e-<'z  } 

A 

rxz  =  'xz  sin(kx~wt) 

bottom: 

with 

?x2  -  -  Qje'52  -  e-1z| 

A 

Tzz  =  Tzz  COS(kX-wt) 

with  _ _ 

Tzz  =  -ksQ  {— X-^y---  e"SZ  +  [6R-2M(X-R)]e"qz  J 

From  Eqs.  15  and  16  we  can  finally  calculate  the  purely  real  energy  flux  in 
the  Scholte  wave. 

The  elastic  energy-flux  density  is  given  by  the  work  done  by  the  stress 
tensor  within  a  unit  time  and  on  a  unit  area  normal  to  the  displacement 
vector: 

e  *  (£.  C,  n)  =  -T s  .  (Eq.  17) 

Accordingly  we  obtain  for  its  non-vanishing  components 
Cw  =  =  ex  cos2(kx  -  wt) 

with 

5  =  ,„k‘M>YHe21"2  .  W|<>Q*Y  _!MML_  e2l"z 

X  S  S  ( X-N ) ( 2X - 1 )2 

water:  (Eq.  18a) 

nw  =  pww  =  ez  sin2(kx  -  wt) 
with: 

e  =  ,K.,k>H*  ^  e2”2  =  wok>  - e2"2 

2  s  2  s  2r1T¥(2X-l)2 
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£  =  -  T  U-T  0) 
S  XX  xz 


=  e'  sin2(kx-o)t)  +  e"  cos2(kx-o)t) 

A  A 


with: 


e;  =  yo3k|Q-Y2(X-R){^T  e~2sz  +  e^z  -  §l|  e"^2  } 


e"  =  yiok2Q2Y  i  8^— -■)  e_2sz  +  [1+2(X-R)] 

x  s  <  (2X-1 )2 


-2qz 


(2X-1 ) 

2AT  /h?  [  1  +2 ( X-R ) ]  „-(s+q)z 
- 2X  -  T  -  e 


} 


bottom: 


(Eq.  18b) 


n  =  -Tzx°  ~  Tzz;  =  ®z  s'*n2(,cx  “  wt) 


with: 


e!  =  rfQ2  /4X/M*  /Mf  e‘2sz  -  (2X-1)2  e"2^z  1 

2  s  2  ( 2X-1 ) 2  {  f 


For  each  point  of  the  interface,  the  term  in  the  braces  {  }  of  the 
amplitude  e^  corresponds  to  the  left-har.d  side  of  the  Scholte  wave  equation 
(Eq.  13)  and  therefore  guarantees  the  continuity  of  the  amplitudes  e^  and 
ez„  As  we  could  further  expect,  the  amplitude  ez  disappears  in  the 
limiting  case  p  =0.  By  taking  the  time  average  of  the  above  results 

W 

(Eq.  18),  the  energy  flux  normal  to  the  interface  vanishes  and  only  the 
component  in  the  propagation  direction  remains: 

A 

-t  ex 

water:  =  £  (Eq.  19a) 

A  A 

e*  +  e" 

bottom:  C  =  ■* ?  x  (Eq.  19b) 

According  to  the  comments  given  in  connection  with  the  displacement 
components  (Eq.  15),  these  mean  energy-flux  densities  also  have  their 
maxima  at  the  interface  and  show  a  very  rapid  frequency-dependent  decay 
with  increasing  distance  from  this  plane,  as  is  shown  in  Fig.  1c. 


For  this  reason,  both  energy  flux  densities  can  usually  be  neglected  for 
distances  |z|  >  A  and  our  assumption  of  treating  the  problem  as  two  half¬ 
spaces  holds  if  both  the  depth  of  the  wate1-  and  that  of  the  sediment  layer 
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are  larger  than  this  limit.  Only  with  an  "extremely- hard"  crystalline  rock 
bottom  may  a  greater  water  depth  be  required  to  guarantee  the  existence  of 
such  a  "pure"  or  "non-modified"  Scholte  wave. 

To  complete  our  mathematical  description  we  write  down  the  expressions  for 
the  total  energy  flux  on  both  sides  of  the  interface  as: 

rw  .  7  C  *  '  ^ 

water:  (Eq.  20a) 

.  „<Q’Y  -MMl - 

s  4.M  (X-N)(2X-1)2 

r  .  Tc  *  dz  ■  nuik2o2Y  /  X(X-RH4X~3)  +  MIM1 

o  5  */H(2X-l)2  VSwf 

bottom:  (Eq •  20b) 

-  (X-R)(4X-1)  +  A-T  >A¥  [1+2(X-R)]  \ 
(/M  +  A^R)(2X-1)  f 


Under  real-world  conditions  the  above-described  wave  type  is  strongly 
modified  of  course  not  only  by  less-pronounced  interfaces  within  the  ocean 
floor  but  also  by  changes  in  parameter  profiles,  inhomogeneities,  and 
material  losses  within  the  resulting  layers  (see  Sects.  1.4  and  2.2). 

Because  of  the  frequency-dependent  penetration  depth  it  is  already  obvious 
that  the  low-frequency  components  of  a  given  Scholte-wave  pulse  are 
affected  more  by  the  bottom  layers,  and  by  the  sea  surface  too  if  the 
occasion  arises,  while  the  high-frequency  components  are  mainly  affected  by 
the  \nhomogeneities  and  the  relative  high  dissipation  of  the  unconsolidated 
top  layers  of  the  sea  floor. 
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1 . 3  Wave-front  Pattern  of  a  Point  Source  Close  to  an  Interface 

Based  on  Cagniard's  rigorous  mathematical  treatment  [ii,39,4cj  (which  does 
not  suffer  from  the  usual  restrictions  required  by  asymptotic  methods)  and 
von  Schmidt's  famous  Schlieren  photos  [41]  (taken  from  model  experiments), 
many  seismic  and  acoustic  papers  have  described  the  different  contributions 
to  the  field  of  an  impulsive  point  source  SQ  located  close  to  a  plane 
interface  separating  two  semi -infinite,  homogeneous,  and  isotropic  elastic 
media.  Unfortunately,  almost  all  these  treatments  are  limited  or  so 
incomplete  that  the  user  who  wants  to  predict  or  interpret  experimental 
results  has  to  make  his  way  through  the  theory  again  and  again.  In  the 
following  we  try  to  describe  the  facts  and  compile  the  consequences  in  a 
clearly  arranged  scheme,  so  that  the  type  and  shape  of  almost  all  possible 
wave  fronts  can  be  quickly  identified.  The  experienced  reader  will  soon 
recognize  that  our  discussion  of  phenomena  makes  use  of  many  more  partial 
contributions  to  the  field  than  those  that  can  be  separated  clearly  by  the 
classical  mathematical  tools  [i?,42]  or  even  more  modern  methods  [22,24,28, 
43  to  51]. 

To  cover  the  most  relevant  cases  we  assume  an  impulsive  point  source  SQ  for 
compressional  (P-)  waves  and  denote  the  two  bulk-  (or  body-)  wave  velocities 
of  the  original  medium  as  c  and  c  and  those  of  the  adjacent  medium  as  c' 

P  b  P 

and  c^.  Besides  this  "hypocentre"  S0>  our  geometrical  considerations  will 
be  based  on  the  fictitious  mirror  source,  S^,  and  on  the  intersection  point 
of  the  connecting  line,  SQS'Q  ,  and  the  interface,  which  is  called  the 
"epicentre"  E,  (nadir  of  SQ).  If  the  expanding  spherical  P-front  impinges 
on  the  flat  interface,  the  line  of  intersection  -  hereafter  called  "trace" 

Tp  -  is  an  enlarging  circle.  This  circle  emanates  from  the  epicentre  E 
with  an  infinitely  high  radial  velocity  Vp  that  decreases  continuously 
until  it  approaches  the  velocity  of  the  incident  spherical  front:  vp  -*•  Cp. 
By  means  of  this  primary  trace,  T  ,  the  incident  front  implements  a  sort  of 
"moving  perturbation"  at  the  interface,  which  has  to  be  balanced  by  an 
elastic  wave  field.  Without  repeating  it  every  time  we  assume  that  all 
outgoing  waves  of  this  resulting  field  fade  away  with  growing  distance  from 
the  source  oQ  ("radiation  condition"). 

To  fulfil  the  transition  conditions  at  the  interface  -  the  continuity  of 
all  displacement  (or  velocity)  and  stress  components  always  requires  the 
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interplay  of  two  compressional  and  two  shear  potentials  -  this  relaxation 
process  is  performed  by  at  least  four  different  radiated  or  guided  waves, 
which  may  be  "pure"  bulk-waves  (body-waves)  or  more  complicated  interface 
phenomena.  The  selection  of  the  wave  types  that  participate  in  the 
transmission  of  the  incoming  energy  is  simply  made  by  the  law  of  causation: 
each  primary  trace  and  all  eventually  emerging  secondary  traces  can  radiate 
or  drag  only  those  fronts  whose  velocity  is  less  than  or  equal  to  the 
momentary  trace  velocity.  Expressed  in  other  words  we  can  say:  the  field 
generated  by  a  source  is  completely  detached  by  radiation  from  its  origin 
or  is  at  least  moving  with  it  if  the  characteristic  propagation  of  this 
field  is  "too  slow"  or  just  "sufficient"  to  follow  the  changes  in  the 
position  of  the  source". 

From  these  fundamentals  we  can  easily  deduce  and  understand  the  well-known 
wave  front  patterns  presented  in  many  papers  without  exhaustive  inter¬ 
pretation:  the  reflected  compressional  front,  PP,  has  always  a  spherical 
shape  -  apparently  it  emanates  from  the  fictitious  mirror  source  - 
while  the  slower,  reflected  shear  front,  PS,  is  the  cap  of  a  hyperbolic 
enclosed  therein.  Strictly  speaking  this  latter  front  is  created  by  "inner 
refraction"  and  not  by  reflection.  As  long  as  the  trace  Tp  is  moving 
faster  than  the  bulk  waves  of  the  second  medium  -  at  least  in  the  early 
stage  of  our  propagation  process  this  condition  is  always  fulfilled  -  the 
refracted  fronts,  PP'  and  PS',  are  also  hyperbolical  caps.  Thus  all 
radiated  wavefronts  intersect  the  interface  at  the  trace  Tp,  with  acute 
angles  av  (correspc  ding  to  the  angles  of  incidence)  that  are  determined  by 
Snells'  law  (see  Fig.  ?a). 

Under  the  condition  c  >  c'  this  quite  simple  set  of  wave  fronts  remains 

r  r 

essentially  unchanged  for  all  time.  During  its  propagation,  all  these 
angles  grow  slowly,  but  only  (the  real  part  of)  a  reaches  its  maximum  tt/2 

r 

in  the  far-field. 

If  the  source  is  located  not  too  far  from  the  interface,  and  additional 
conditions  concerning  the  elastic  parameters  are  fulfilled,  free  interface 
waves  form  the  "tail"  of  the  resulting  elastic  wave  field.  As  mentioned 
above,  the  excitation  of  these  slowest,  guided  contributions  is  due  to  the 
inhomogeneous  character  of  the  incident  spherical  front,  which  can  be 
demonstrated  by  the  classical  mathematical  treatment  of  boundary-value 
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problems  in  homogeneous  and  isotropic  media.  The  first  step  towards  the 
solution  usually  consists  in  adapting  the  original  representation  of  the 
wave  field  to  the  geometry  of  the  boundary  or  interface.  This  modification 
has  to  be  done  in  such  a  way  that  the  resulting  separation  of  variables 
allows  a  suitable  fulfilment  of  the  transition  conditions.  In  our  quite 
simple  case  of  a  flat  interface  these  adapted  representations  are  given  by 
the  already-mentioned  contour  integrals  derived  by  Sommerfeld  and  Weyl 
[3  to  9].  These  expressions  decompose  an  incoming  spherical  or  cylindrical 
field  into  an  infinite  sum  of  cylindrical  or  plane  waves  with  complex  wave 
vectors  or  complex  angles  of  incidence.  The  imaginary  part  of  these 
vectors  or  angles  describes  the  inhomogeneity  of  the  individual  element  in 
relation  to  its  distance  from  the  interface.  Because  they  represent 
crosswise  evanescent  waves  of  a  subordinate  geometry,  these  components  are 
provided  with  features  that  are  similar  to  the  free  interface  waves  to  be 
excited.  In  our  special  case  the  genuine  type  is  characterized  by 
cylindrical  phase  surfaces  orthogonal  to  the  "guiding"  interface  and  plane 
amplitude  surfaces  parallel  to  it  and  thus  by  a  restricted  penetration 
depth  inside  both  media.  Synthesizing  the  incoming  field  from  those 
inhomogeneous  elements  we  get  integrals  that  resemble  the  well-known 
Fourier  or  Fourier-Bessel  transformation  [52,53].  From  the  subsequent 
integration  along  a  suitable  path  on  the  admissible  Riemann  sheet,  we  can 
work  out  in  general  the  contributions  of  one  saddle-point  (reflected  wave), 
up  to  four  branch-points  (refracted  bulk  waves),  and  up  to  two  major  poles 
(free  interface  waves)  [11,17,23,42,43,  46  to  51], 

The  facis  so  far  outlined  characterize  completely  the  simplest  case  A  of 
Taole  2,  assuming  cp  >  cp.  In  this  scheme  we  have  classified  all  outgoing 
wave  fronts  as  "non-homogenous"  because  of  the  complicated  amplitude 
weighting  of  the  reflected  and  refracted  fronts.  Their  deviation  from 
homogeneity  is  easy  to  understand  if  one  bears  in  mind  that  the  reflection 
and  transmission  coefficients  depend  on  the  angle  of  incidence  according  to 
the  way  that  the  ensemble  of  these  outgoing  waves  forms  the  "radiation 
pattern"  of  the  insonified  interface  region. 

To  simplify  Table  2  we  have  characterized  all  wave  fronts  by  only  the  last 
and  most  significant  letter  of  their  usual  seismic  classification.  Under 
the  condition  cp  ■>  cp  >  c^  of  case  B  in  Table  2  there  occurs  a  moment  at 
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which  the  decreasing  trace  velocity  Vp  reaches  and  passes  the  value  of  the 
congressional  velocity  Cp  in  the  second  medium.  At  this  intermediate  stage 
of  the  resulting  field  the  real  part  of  angle  cip  of  the  refracted  front 
PP1  adopts  its  (real)  maximum  tt/2  (see  Fig.  2b)  and  continues  with  complex 
values.  Thus  the  refracted  front  PP',  having  a  hyperbolical  shape  until 
now,  is  converted  into  a  half-ellipsoid  that  ends  orthogonally  at  the 
interface  and  expands  with  its  characteristic  velocity  c^.  Accordingly, 
this  modified  front  PP1  has  to  be  disconnected  from  the  primary  trace  Tp 
(indicated  by  "o"),  which  is  dropping  back  continuously.  In  term  of  ray 
theory  this  moment  corresponds  to  the  reaching  and  exceeding  of  the 
critical  angle  of  incidence  ap  cr- t  with  the  "disappearance"  of  the 
refracted  front  PP1  and  the  corresponding  intensification  of  the  reflected 
front  PP.  To  fulfil  the  later  transition  conditions,  the  field  emanating 
from  the  relative  "slow"  trace  Tp  has  to  be  "completed"  again  by  the  aid  of 
a  new  compressional  component  in  the  second  medium.  Because  it  is  dragged 
but  not  radiated  by  the  (primary)  trace  T  ,  this  attached  "near-field"  is 
created  by  a  guided,  inhomogeneous  wave  F*P'  belonging  to  a  class  of 
evanescent  fields  that  are  frequently  called  "flank-waves"  (see  Fig.  2c). 

In  Table  2  these  coupled  annular  flank  waves  are  indicated  by  downward 
pointing,  dotted  arrows.  Due  to  their  low  energy  level  they  have  almost  no 
importance  for  seismic  investigations,  but  together  with  the  free  interface 
waves  they  exercise  an  influence  on  the  fine-structure  of  totally-reflected 
beams  [57,  60] ,  which  becomes  quite  complicated  close  to  the  so-called 
Rayleigh  angle  [61].  On  the  other  hand,  the  relative  "fast"  front  PP1  has 
to  create  its  own  (secondary)  trace  T^  (indicated  by  "0")  with  at  least 
three  additional  contributions  to  the  field.  For  historical  reasons  those 
secondary  bulk-wave  traces,  such  as  T^,  are  designated  somewhat  misleadingly 
as  "lateral  waves".  As  it  propagates  with  a  higher  velocity  than  all  other 
wave  types,  trace  can  easily  satisfy  the  transition  conditions  by 
radiating  residual  "pure"  bulk  waves.  Because  they  are  dragged  by  and 
tangent  to  the  respective  fronts  of  the  same  type  (emanating  from  T  )  these 
new  wave  fronts,  PP'P,  PP'S,  ai  PP'S',  have  a  conical  shape  (see  Fig.  2c). 
They  are  the  famous  "head  waves"  that  were  long  an  unsolved  problem  in 
seismology  until  Cagniard  succeeded  in  explaining  them  theoretically  and 
von  Schmitt  proved  their  existence  experimentally.  In  our  scheme  they  are 
designated  in  the  last  column  by  upward  pointing  arrows  emanating  from  the 
corresponding  secondary  traces.  As  these  head  waves  continuously  draw 
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their  energy  from  the  pertinent  lateral  waves  (traces),  it  is  obvious  that 
their  amplitudes  must  always  be  greater  in  sections  previously  radiated 
than  in  subsequent  ones.  This  fact  is  indicated  schematically  in  Fig.  2c 
by  the  line  thickness  of  the  three  corresponding  fronts.  Thus  their 
amplitude/distance  characteristic  is  apparently  a  rising  function  with 
increasing  distance  from  the  interface  and  reaches  its  maximum  at  the 
junction  with  the  respective  bulk-wave  front.  Such  a  head  wave  provokes 
the  one-sided  radiation  loss  of  the  pseudo-interface  waves  described  in  the 
previous  Section. 

Under  the  condition  Cp  >  c^  >  cp  °f  case  C  in  Table  2  the  refracted  shear- 
front  PS1  is  also  converted  into  a  half-ellipsoid  and  thus  creates  its  own 
(secondary)  trace  T^  subsequent  to  the  separation  of  Tp  and  T^.  This 
process  corresponds  to  the  reaching  and  exceeding  of  a  second  critical 
angle  of  incidence  a$  ^  characterized  by  the  "disappearance"  of  the 
ultimate  refracted  front  PS'  and  thus,  in  fact,  by  the  "total"  reflection 
of  the  incoming  front  P.  Due  to  its  intermediate  velocity  this  third  trace 
T^  drags  two  head  waves,  PS'P  and  PS'S,  in  the  first  medium  but  only  a 
compressional  flank  wave  PS'P' in  the  second  medium. 

After  this  exhaustive  discussion  we  can  easily  translate  the  most  compli¬ 
cated  case  C  of  our  scheme  into  the  corresponding  wave-front  pattern  of 
Fig.  3,  where  it  is  seen  to  be  characterized  by  the  occurrence  of  five  head 
waves  and  three  flank  waves.  In  general,  the  interesting  free  interface 
waves  cannot  be  incomorated  in  such  a  clear  illustration  with  the  same 
surety,  because  their  number  (0,1,  or  2)  and  type  has  to  be  investigated 
for  each  individual  case.  The  shape  of  their  phase  surface  is  cylindrical 
if  they  form  the  "slowest"  contribution  to  the  field  (genuine  interface 
waves)  otherwise  composed  from  an  upen  cone  in  the  "fast"  medium  and  an 
almost  closed  head-wave  cone  in  the  adjacent  "slow"  medium  (pseudo¬ 
interface  waves)  [46] . 

It  should  be  mentioned  that  the  above  description  of  the  different  wave 
fronts  is  restricted  to  the  most  important  phenomena.  Cagniard  for 
instance  has  shown  that  all  hyperbolical  and  elliptical  fronts  are  followed 
by  blurred  spherical  fronts  of  the  same  type,  which  are  like  a  "shadow"  of 
their  early  stage.  These  weak  "archetypes"  create  additional  obliterated 
traces  at  the  interface  with  the  corresponding  remainders  of  head  and  flank 
waves.  Those  very  weak  flank  waves  are  called  "surface  waves  of  the  second 
type-1  by  Cagniard,  but  a  much  less  confusing  name  for  them  would  simply  be 
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FIG.  3  THE  MOST  COMPLICATED  WAVE-FRONT  PATTERN  AT  THE  INTERFACE 
OF  TWO  SCUDS  (see  Case  C  of  Table  2) 
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"smeared  flank  waves".  In  his  study  of  leaking  interface  waves,  Phinney 
[32]  describes  a  so-called  pseudo-P  wave  that  may  appear  as  an  appendage  to 
the  front  of  the  P-wave  or  may  become  a  distinct  signal  between  the  arrival 
of  the  P  and  S  fronts.  Its  particular  velocity  always  adopts  that  value 
that  minimizes  the  losses  by  conversion  to  radiated  head  waves. 

At  first  glance,  the  above-described  occurrence  of  several  non-homogeneous 
and  inhomogeneous  wave  fronts  with  a  different  shape  may  remind  the  reader 
of  apparently  similar  phenomena  in  crystal  optics,  such  as  double  and 
conical  refraction,  for  example.  A  cursory  examination  reveals,  however, 
the  complete  difference  in  the  underlying  physics.  With  the  matter- 
connected  elastic  waves  this  division  into  several  wave  fronts  is  already 
the  consequence  of  all  possible  independent  oscillations  of  the  microscopic 
elements  relative  to  the  regular  lattice  or  the  molecular  cluster  of  the 
medium  and  relative  to  the  shaping  interfaces  of  the  body.  With  the 
matter-independent  electromagnetic  waves  consisting  merely  of  one  basic 
type  with  transverse  field  vectors  {at  least  when  far  from  the  source)  this 
splitting  effect  is  due  to  the  relative  orientation  of  the  radiated  field 
vectors  and  the  anisotropic  microstructure  of  the  exposed  medium  (polariza¬ 
tion  effects) . 


1.4  Sound  Propagation  in  a  Stratified  Oceanic  Lithosphere 
and  a  Wedge-shaped  Water  Layer  on  top 

By  discussing  the  propagation  of  elastic  waves  inside  unlimited,  homogeneous, 
and  isotropic  media,  and  along  the  plane  interfaces  that  separate  them,  we 
have  covered  the  fundamentals  of  more  complicated  phenomena  occurring  in 
finite  bodies  composed  of  different  media.  If  the  sound  field  interacts 
with  more  than  one  of  the  interfaces  that  bounds  the  medium  the  totality  of 
all  internally  reflected  and  refracted  components  may  create  other  mixed- 
wave  types  that  are  tuned  to  the  corresponding  body  as  a  whole  and 
synthesize  its  natural  vibrations  or  modes.  As  this  tuning  is  based  on 
constructive  interference  in  the  resulting  field,  the  dimensions  of  the 
body  have  to  be  large  enough  with  respect  to  the  sound  frequencies  above 
certain  cut-off  limits.  Thus  a  free  plate  represents  a  wave  guide  for 
(relatively  high-frequency)  Rayleigh  waves  and  an  infinite  number  of  normal 
modes;  this  is  usually  given  the  collective  name  of  Lamb  waves  [ 6,17,21 ,  62 
to  6b],  It  should  here  be  mentioned  that  Lamb  was  the  first  to  tackle 
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successfully  the  half-space  problem  in  a  classical  paper  [16]  that  already 
included  a  very  useful  mathematical  description  of  Rayleigh-wave  pulses. 

If  these  mixed  modes  in  plates  provoke  symmetrical  deformations  with 
respect  to  the  middle  plane  (due  to  transverse  contraction)  they  are 
called,  by  analogy,  "longitudinal"  or,  more  correctly,  "di latational "  or 
"extensional "  waves  (sv),  while  their  antisymmetric  counterparts  (being 
essentially  shear  phenomena)  form  the  well-known  family  of  "bending"  waves 

<*»>• 

In  general,  the  numbering  of  those  modes  reflects  the  number  of  zero- 
crossings  passed  through  by  the  amplitude  distribution  along  the  cross- 
section  of  the  wave  guide.  For  both  symmetries  the  zeroth-order  mode 
occupies  a  special  position:  it  adopts  finite  values  of  phase  velocity  in 
the  low-frequency  range  (no  cut-off)  and  the  high-frequency  limit  of  its 
dispersion  curve  is  determined  by  the  speed  of  the  Rayleigh  wave,  as  we 
have  plotted  qualitatively  in  Fig.  4. 


FIG.  4  DISPERSION  CURVES  FOR  THE  SYMMETRICAL  (sv)  AND  ANTISYMMETRIC  (av> 
PROPAGATION  MODES  IN  A  FREE  PLATE 
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Because  it  represents  a  sort  of  elastic  counterpart  of  the  electromagnetic 
skin  effect,  the  latter  feature  has  lead  to  many  applications  in  ultra¬ 
sonics  using  the  Rayleigh  wave  as  a  tool  for  non-destructive  testing  (see 
e.g.  [66,67])  and  signal  processing  (see  e.g.  [68,69]).  By  increasing  the 
frequency  (or  the  plate  thickness)  the  dispersion  curves  of  all  higher- 
order  modes  start  with  an  infinite  phase  velocity  at  the  corresponding  cut¬ 
off  frequency  and  drop  to  the  shear  speed  of  the  plate  medium. 

A  liquid  layer  on  a  solid  sub-stratum  can  support  only  a  family  of  pure 
compressional  modes  that  behave  similarly  to  the  extensional  waves  of  the 
plate,  but  cause,  in  general,  asymmetric  deflections  of  the  layer  [5,70] . 

The  zeroth  mode,  also,  has  no  cut-off  frequency  and  its  phase  velocity  is 
identical  with  the  Rayleigh  speed  in  the  range  of  extreme  low  frequencies 
(or  for  very  thin  layers)  and  approaches  the  Scholte  speed  for  extreme  high 
frequencies  (or  very  thick  layers)  [71].  Under  the  quite  rare  condition 
that  the  sound  speed  of  the  liquid  exceeds  both  the  bulk-wave  velocities  of 
the  solid  substratum,  no  other  modes  of  higher  order  can  be  propagated 
through  the  liquid  layer.  Usually  the  sound  speed  of  the  liquid  is  smaller 
than  one  or  both  bulk  velocities  of  the  underlaying  solid,  so  that  an 
infinite  number  of  higher-order  modes  can  be  excited,  as  indicated 
qualitatively  in  Fig.  5.  Going  up  the  frequency  scale  (or  increasing  the 
layer  thickness)  the  dispersion  curves  of  all  these  higher-order  modes 


C  PHASE  VELOCITY 


FIG.  5  PHASE-VELOCITIES  OF  THE  LOWEST  MODES  IN  SHALLOW  WATER 
OVER  A  LIQUID-LINE  OR  ROCK-LIKE  BOTTOM 
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start  with  the  bulk-velocity  of  the  solid  (at  the  corresponding  cut-off 
frequencies  f  )  and  drop  to  the  sound  speed  cQ  of  the  liquid.  The  bulk 
velocity  cv  is  identical  with  the  congressional  speed  cp  under  the 
condition  (A):  cp  >  cw  »  c$  (such  as  shallow  water  overlying  a  "very 
soft"  sedimentary  sea  floor),  or  identical  with  the  shear  speed  cg  under 
the  condition  (B):  cp  >  cg  >  cw  (such  as  s  allow  water  overlying  "very 
hard"  rock-bottom).  These  two  possibilities  can  be  easily  understood  in 
terms  of  mode  synthesis  by  constructive  interference. 

In  Fig.  6  we  have  therefore  plotted  qualitatively  the  absolute  value  of  the 
reflection  coefficient  for  both  conditions  A  and  B  against  the  grazing 
angle  Bw  =  7  -  a  ,  i.e.  the  complement  of  the  angle  of  incidence  a. 

[5,72].  The  smooth  dotted  curve  takes  into  account  material  losses  and 
therefore  approaches  the  familiar  experimental  results.  While  condition  A 
offers  the  possibility  of  constructive  interference  for  all  angles  3 

W 

remaining  below  3p  cr^t»  condition  B  can  do  this  only  below  the  usually 

smaller  angle  3C  and,  in  principle,  for  a  small  region  around  & 

5  ,cn  L  p  ,cn  L 

The  narrow  width  of  the  peak  at  Bp  cr^t  and  its  suppression  under  realistic 
conditions  indicates  that  this  interval  is  of  no  importance  in  practice. 

Thus  the  upper  velocity  limit  of  the  higher-order  modes  is  determined  by  cp 
for  condition  A  and  by  cs  for  condition  B. 

So  far,  the  "mode  picture"  has  served  as  a  very  familiar  and  convenient 

scheme  to  explain  and  to  illustrate  the  consequences  of  our  shift  to  the 

low-frequency  regime.  At  this  point,  however,  we  should  draw  the  reader's 

attention  to  the  fundamentals  of  the  normal-mode  theory.  Only  condition  B 

deals  with  a  real  eigenvalue  problem  and  thus  justifies  our  reasoning  for 

those  discrete  modes  and  well-defined  cut-off  conditions.  Condition  A 

already  belongs  to  the  more  general  complex  eigenvalue  problems  that  do  not 

permit  strict  mode  assignments  except  for  the  interesting  Scholte  wave, 

which  always  represents  a  discrete  propagation  mode,  as  we  have  shown  in 

Sect.  1.2  (zeroth  mode  with  respect  to  the  above-mentioned  classification 

of  the  water  wave).  With  our  limitation  to  "very  soft"  sediments  (c  »  c  ) 

w  s 

we  presumed  tacitly  a  quasi-discrete  eigenvalue  problem  corresponding  to  an 
almost  liquid-like  bottom.  As  usual,  the  most  common  and  realistic 
environmental  conditions  lie  between  (cw  »  c$;  "normal"  bottom-types,  such 
as  "relatively  hard"  sediments  or  "relatively  soft"  sedimentary  rocks)  and 
therefore  create  a  more  complicated  and  much  less  transparent  field 


FIG.  6  TYPICAL  REFLECTION  COEFFICIENTS  AT  A  VER 
SEA  FLOOR  (A)  AND  AT  A  VERY  HARD  ROCK-BO 
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composition.  In  Fig.  7  the  dotted  curves  are  intended  to  demonstrate  that 
the  pertinent  blurred  interference  patterns  still  behave  like  a  non- 
systematic  and  less  pronounced  mode  structure.  We  will  therefore  continue 
to  use  the  vivid  mode  picture  and  use  arguments  based  on  those  discrete 
phenomena  to  describe  qualitatively  the  underlying  physics. 

The  two  simple  examples  above  should  be  the  starting  point  from  which  to 
obtain  at  least  some  impression  of  the  quite  complicated  wave  propagation 
in  the  earth's  crust.  The  ocean-floor  environment  that  is  of  particular 
interest  to  is  in  this  report  is  an  irregularly-layered  structure  of 
inhomogeneous,  anisotropic,  and  dissipative  media  with  pronounced  parameter 
profiles.  For  example,  densities  and  bulk  wave  velocities  increase  with 
depth  and  losses  decrease  simultaneously  (see  Sect.  2.2).  All  internal 
interfaces  impose  additional  transition  conditions  on  the  elastic  field  and 
they  may  support  different  types  of  interface  waves  (see  Sect.  1.1  and  [73]). 
Furthermore,  some  of  the  layers  may  be  distinct  wave  guides  for  certain 
trapped  modes  [74/75],  such  as  Love  waves  in  sedimentary  layers  [5,76,77]. 

While  the  general  influence  of  bottom  stratification  on  the  reflection 

coefficient  was  carefully  explored  by  underwater  acousticians  and 

geophysicists  in  the  past  (see  e.g.  [78  to  82]),  it  was  only  a  short  time 

ago  that  Hawker  [83,84]  demonstrated  how  the  specific  influence  of 

interface  waves  resulted  in  a  pronounced  minimum  of  this  coefficient 

occurring  next  to  grazing  incidence,  i.e.  for  small  angles  6  .  Due  to  the 

w 

greater  variety  of  seismic-wave  types,  the  dispersion  curve  of  the  "zeroth 
mode"  of  the  system  may  be  split  into  different  branches  forming  several 
"sub-modes"  characterized  by  the  pronounced  interfaces.  The  high-frequency 
limit  of  each  of  these  separated  branches  should  be  strongly  related  to  the 
velocity  of  a  Scholte  or  Stoneley  wave  at  the  upper  margin  of  a  pronounced 
top  layer,  while  the  corresponding  low-frequency  limit  should  be  linked  to 
a  somewhat  higher  velocity  characterizing  a  sort  of  "modified  Rayleigh 
wave"  of  the  underlaying  stratified  half-space.  Generalizing  this  dynamic 
elasticity  of  simple  structures,  the  phase  velocities  of  higher-order  modes 
should  vary  between  the  most  significant  bulk-wave  velocities  of  the  stack. 
Thus  we  can  expect  that  their  upper  limit  at  the  low-frequency  cut-off  is 
determined  by  the  sub-bottom,  while  their  lower  limit  in  the  high-frequency 
range  is  a  sort  of  "mean  value"  that  refers  to  the  top  layers.  Obviously 
the  exitation  of  those  resulting  modes  that  insonify  the  complete 
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FIG.  7  PHASE-VELOCITIES  OF  THE  LOWEST  MODES  IN  SHALLOW  WATER  OVER 
A  RELATIVE  "SOFT"  ROCK-BOTTOM 
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structure  and  of  those  trapped  modes  that  are  restricted  to  certain  layers 
or  interfaces  depends  on  the  given  frequency  scale  and  last,  but  not  least, 
on  the  type  and  position  of  the  source.  The  successful  application  of  the 
normal -mode  model  for  many  sonar  purposes  [5,  85  to  90]  shows  that  very 
simple  bottom  models  may  already  be  sufficient  to  make  allowance  for  the 
influence  of  the  sea  floor  on  water-borne  sound  in  the  medium-  and  high- 
frequency  regime.  However,  in  the  range  of  low  and  very  low  frequencies  we 
have  to  proceed  to  seismic  models  of  the  earth's  crust,  because  the  water 
depth  and  the  thickness  of  many  bottom  layers  may  be  comparable  to,  or  even 
smaller  than,  one  wavelength  of  the  transmitted  sound  field.  A  reasonable 
model  for  such  infrasonic  propagation  problems  is  already  given  by  a 
horizontally-stratified  half-space  composed  of  homogeneous  and  isotropic 
layers  characterized  by  constant  parameters.  A  progressive  refinement  of 
the  idealized  layering  provides  us  moreover  with  the  possibility  of 
approximately  shaping  the  typical  parameter  profiles  if  needed  (see  Sect. 
1.3). 

Other  very  common  features  of  sound  propagation  in  shallow  water  are  the 
coupling  or  conversion  effects  of  different  modes  belonging  to  the  same 
wave  type  or  even  to  quite  dissimilar  propagation  phenomena.  Those 
interactions  may  already  happen  in  regularly-layered  systems  with  range- 
dependent  parameter  profiles  (see  e.g.  [26, 70]),  but  they  are  typical  for 
structures  with  varying  layer  thicknesses  (see  e.g.  [91  to  93])  and/or 
range-dependent  parameter  profiles  (see  e.g.  [94]). 

To  end  this  discussion  of  the  influence  of  layering  on  wave  propagation  we 
would  therefore  like  to  drop  a  hint  on  the  effects  of  a  sloping  bottom.  As 
they  propagate  towards  the  apex  of  the  wedge-shaped  liquid  layer,  the 
higher-order  modes  of  the  "water  wave"  are  "stripped"  one  after  another. 

On  each  occasion  the  energy  is  partly  transferred  to  the  next  lower-order 
mode  (mode  conversion),  partly  reflected  back  into  deeper  water,  and  partly 
refracted  into  the  substratum  in  the  form  of  a  well-defined  beam  consisting 
mainly  of  radiated  P-waves  ("soft"  bottom)  or  SV-waves  ("hard"  bottom). 
While  these  coupling  effects  (conversion  from  multi-mode  to  single-mode 
propagation)  and  radiation  characteristics  of  tapered  wave  guides  have 
already  been  studied  carefully  in  connection  with  guided  micro-waves  and 
light  waves  (see  e.g.  [95-97]) , their  counterparts  in  underwater  acoustics 
have  attracted  little  attention  in  the  past.  With  the  growing  interest  in 
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seismic  sensing  the  latter  problems  can  now  be  expected  to  experience  a 
revival.  Based  on  the  results  with  the  less-complicated  electromagnetic 
fields,  Coppens  and  Sanders  [98]  and  Odom  et  al  [99]  started  to  model  theo¬ 
retically  and  experimentally  the  shallow-water  edge  with  different,  quite- 
small,  taper  -  or  "beach"  -  angles.  In  accc  dance  both  with  simplified 
theories  and  model  experiments,  they  demonstrated  that  the  acoustic  energy 
is  rapidly  dumped  into  the  bottom  before  the  first  mode  reaches  the 
pertinent  cut-off  depth.  The  refracted  congressional  or  shear  wave  in  the 
substratum  always  formed  a  well-defined  beam  (narrower  with  decreasing 
slope  or  taper  angle)  but  their  predicted  depression  angle  (about  15°  for 
most  realistic  parameter  combinations)  always  exceeded  the  measured  result 
(below  10°).  Modelling  this  problem  with  a  more  sophisticated  mathematical 
tool  (Parabolic  Equation)  Jensen  and  Kuperman  [ioo]  reached  better 
agreements  with  those  experimental  data  (for  a  liquid-like  bottom)  and  also 
showed  the  corresponding  refraction  phenomena  for  some  of  the  higher-order 
modes . 

The  top  layers  of  the  ocean  bed  mostly  consist  of  unconsolidated  sediments 
(condition  A),  which  is  why  the  primary  conversion  or  bottom  coupling  of 
the  suppressed  modes  usually  creates  a  P-wave  in  the  sea-floor.  The 
pertinent  beam  then  suffers  repeated  and  more  complicated  reflection-  and 
refraction-splittings  due  to  deeper  interfaces,  and  simultaneously 
undergoes  a  strong  broadening  as  a  consequence  of  the  different  upward- 
refracting  parameter  profiles  and  inhomogeneities  within  the  layers  (see 
Sect.  2.2).  This  latter  comment  on  the  stripping  of  all  higher-order  modes 
underlines  again  the  special  position  of  the  Scholte  wave  or  its  modifica¬ 
tions  (see  Sect.  3.2):  that  they  are  not  only  free  from  those  radiation 
losses  but  may  even  profit  from  them  (see  e.g.  83,84]). 


1 . 5  Some  Analytic  Approaches  and  Model  Experiments  Regarding 
Annular  Scholte-wave  Pulses 


As  mentioned  in  Sect.  1.1,  in  the  1930s  Cagniard  [ii]  solved  the  problem  of 
an  impulsive  point  source  located  close  to  a  plane,  solid/solid  or 
liquid/solid  interface.  His  famous  "method"  of  deriving  the  closed-form 
solutions  for  all  ranges  (see  also  [39,40])  is  essentially  based  on  Lamb's 
work  [17]  and  was  developed  independently  and  almost  simultaneously  by 
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Smirnow  and  Sobolev  [101,102].  However,  due  to  the  sophistication  of  the 
mathematics  involved,  these  results  never  became  a  practical  tool  to  solve 
geophysical  or  acoustical  problems. 

Thus  after  World-War  II  a  few  physicists  tried  to  find  more  practicable 
approaches  and,  especially  in  the  years  1960-65,  some  efforts  were  devoted 
to  the  derivation  of  much  simpler  mathematical  descriptions  and  to  their 
"testing"  with  the  aid  of  model  experiments  in  the  laboratory  (e.g.  [103]). 
It  is  obvious  that  the  results  of  those  small-scale  experiments  can  merely 
reflect  some  very  pronounced  features  of  real-world  data  but  they  may 
provide  at  least  a  very  helpful  cross-check  on  the  theoretical  set-up  and 
give  some  confidence  in  applying  it  under  more  realistic  conditions. 

Before  discussing  some  of  these  efforts  in  detail  we  would  like  to  review 
the  pertinent  mathematics  in  general.  The  wave  propagation  of  interest 
along  an  interface  is  characterized  by  up  to  four  branch  points  of  the 
pertinent  integrands,  corresponding  to  a  maximum  of  four  "undisturbed"  bulk 
waves.  To  restore  the  uniqueness  of  the  complex  plane  of  integration  we 
have  therefore  to  choose  up  to  four  branch  lines  or  "cuts"  that  divide  the 
ambiguous  Riemann  surface  of  the  integrals  into  a  maximum  of  five  "sheets" 
or  "leaves".  Due  to  the  different  ways  of  making  these  cuts,  and  to  the 
resulting  definition  of  those  sheets,  one  has  often  several  possibilities 
for  calculating  the  different  contributions  to  the  integrals.  The  final 
choice  is  usually  determined  by  the  requirement  for  a  fast  convex nence  and 
by  a  need  to  facilitate  the  evaluation.  For  a  fixed  source/receiver 
configuration  those  joined  sheets  can  be  interpret’d  quite  plausibly  as  the 
subsequent  "time  windows"  of  the  synthetic  seismogram  (being  separated  by 
the  arrivals  of  the  bulk  waves)  and  thus  their  contributions  can  be 
interpreted  as  the  pertinent  sections  of  the  calculated  signal.  The  only 
sheet  without  any  contribution  and  thus  "without  physical  meaning"  is 
usually  classified  as  the  "lowest"  one  corresponding  to  the  period  between 
the  onset  of  the  source  and  the  first  arrival  at  the  receiver  position. 

The  genuine  interface  waves  of  interest  are  characterized  by  poles  on  the 
"top  sheet"  that  contribute  to  the  "tail"  of  the  signal,  while  the  pseudo¬ 
interface  waves  result  from  poles  on  "lower  sheets"  between  the  mentioned 
extremes. 
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Despite  the  fact  that  the  integration  has  to  be  performed  on  the  top  sheet, 
those  lower  sheet  poles  may  have  an  influence  on  the  synthetic  seismograms 
by  their  proximity  to  the  cuts.  With  respect  to  the  more  practicable 
approaches  that  reveal  this  correlation  between  the  mathematical  and  the 
physical  aspects,  the  three-part  paper  by  Roever,  Vining  and  Strick  [24] 
and  the  publication  by  Phinney  [32]  may  be  called  milestones  in  the  field 
of  elastic  interface  waves. 

To  simplify  the  mathematics,  Strick  reduced  the  problem  to  two  dimensions, 

assuming  a  delta-function-excited,  horizontal,  line  source  instead  of  a 

similar  point  source.  As  the  resulting  field  terms  of  these  sources  differ 

at  large  ranges  by  essentially  an  amplitude  factor  and  a  half-order 

differentiation,  such  a  line  source  may  be  regarded  as  a  point  source  with 

an  excitation  between  a  delta-  and  a  step-function.  The  chosen  line  source 

is  therefore  a  good  simulation  of  a  detonating  chemical  charge  or  of  an 

electrical -spark  discharge,  which  are  the  commonest  source  devices  for 

field  studies  and  laboratory  tests.  The  only  substantial  adulteration 

results  from  the  characteristic  "tail -effect"  of  a  line  source,  which  can 

be  easily  understood  from  Huyghen's  principle.  (The  source  cannot  be 

blanketed  off  and  thus  cannot  be  weighted  by  an  amplitude  assignment.) 

Assuming  two  adjacent  half-spaces,  Strick  calculated  the  pressure  response 

and  the  particle  velocities  at  a  water/solid  interface  corresponding  to  a 

"relatively  soft"  (condition  I:  cc  <  c..)  and  an  "extremely  hard" 

s  w 

(condition  II:  c  >  c.J  sea  floor.  In  both  cases  he  proved  the  excitation 
of  a  relatively  strong  Scholte-wave  pulse,  and  demonstrated  impressively 
that  the  already-mentioned  suporession  of  the  pseudo-Rayleigh  wave  at  a 
soft  ocean  bottom  (for  the  physical  interpretation  see  Sect.  1.2)  is  due  to 
the  off-axis  position  of  the  pertinent  pole.  In  addition,  Strick  has 
demonstrated  that  the  real  part  of  this  pseudo-Rayleigh  pole  is  practically 
identical  with  the  true,  real -Rayleigh  pole  and  thus  with  the  zero- 
crossings  of  the  pertinent  signals  propagating  with  the  same  velocity. 

Figure  8  compares  two  examples  from  Strick's  theoretical  results  with  the 
pressure  records  of  a  transient  point  source  (spark  gap)  made  by  Roever  and 
Vining.  The  agreement  is  indeed  outstanding,  and  as  was  to  be  expected 
(from  our  theoretical  findings  in  Sects.  1.1  and  1.2)  the  Scholte-wave 
pulse  is  more  compressed  under  "very  hard"  bottom  conditions  but  also  more 
difficult  to  separate  from  the  direct  P-wave  arrival  in  the  water  column. 
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FIG.  8  COMPARISON  OF  STRICK'S  MODELLING  RESULTS  WITH  THE  LABORATORY 
MEASUREMENTS  BY  ROEVER  AND  VINING  [23] 
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Stride's  mathematical  treatment  of  a  given  source/interface/receiver 
constellation  already  demonstrates  the  decisive  role  of  the  so-called 
"period-function"  F  and  "response-function"  6,  which  are  given  by: 


F  =  /jmST  (Sr  -  sL) 


(Eq.  21a) 


G  =  /FN  (Sr  +  SL)  , 


(Eq.  21b) 


with:  X,  N  defined  by  Eq.  11. 

Sp,SL  right-and  left-hand  side  of  Eq.  13. 


For  an  exponentially-decaying  source  pulse  (pressure  p  or  displacement 
u,w),  such  as:  t 

X0  =  30  e  ;  (Eq.  22) 


with  B  .  1  =  constant  real  factors, 
o  o 

Phinney  has  shown  that  for  T  ^  0  and  BQ  -►  °°  (i.e.  XQ  -*■  6(t))  the  excited 
wavelets  of  the  genuine  or  pseudo-interface  waves  can  be  described  by: 


BiB2(u)B3(r,t) 

(„  m  /  a  A„  m - cos  fL(u>,r, t) 

^n,m  /  ~  n,m  r—?  nv  ' 

/7=0  v  r 


(Eq.  23) 


P)  (0 J  (0 

u  >  for  n  =  <  1  >  and  m  =  <  1 
w  1  2 


An  m  =  complex  amplitude  factor 


Bi  =  lim  (BQT0)  constant  source  parameter 
o*  ^ 


C(F,G)  =  i  (  — —  =  B2(to)  eili;(to) 

9F/3o)  /(i)=kc 


system's  "excitation  function". 

D(r,t)  =  (r-cw/f  t)  -  ih  A-  ^  =  B3(r,t)  e1<^r,t) 


signal's  "pulse-envelope". 
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ftn(<i>,r,t)  =  Df>(to)  +  n  ^  4> ( ^ t ) 
signal's  "phase  angle". 

The  most  prominent  features  of  the  general  expression  of  Eq.  23  are 
determined  by  the  expected  term  for  cylindrical  spreading,  by  a  charac¬ 
teristic  weighting  function  C  of  the  system  (having  the  well-known 
features  of  a  reflection  or  transmission  coefficient),  and  by  a  sort  of 
"carrier"  having  the  variable  frequency  and  being  shaped  by  the 
"amplitude-modulation"  D.  As  the  results  for  u  and  w  are  very  similar  but 
without  any  point  of  symmetry  (which  occurs  in  such  two-dimensional 
problems  as  the  above-discussed  approach  by  Strick)  and  are  phase-shifted 
against  each  other,  the  resulting  particle  orbit,  and  thus  the  hodographs 
of  particle  velocity  and  acceleration,  are  generalized  lemniscates  (see 
also  [11]). 

In  agreement  with  our  monochromatic  findings  in  Sect.  1.2,  the  broadband 
results  of  Strick  and  Phinney  demonstrate  that  the  effect  of  removing 
source  or  receiver  from  the  interface  is  to  suppress  the  higher  frequencies 
of  the  excited  interface  wave.  If  the  radiated  pulse  shape  and  the 
receiver  depth  are  well  known,  the  spectrum  of  the  detected  signal  reflects 
the  distance  of  tne  source  from  the  interface. 

Under  real-world  conditions  this  feature  is  certainly  of  secondary 
importance  because  of  the  effects  of  dispersion  (see  Sect.  1.4)  and 
absorption  (see  Sect.  2.2). 

Between  1963  and  1965,  Spitznogle,  McLeroy,  and  co-workers  extended 
Strick's  theoretical  work  and  the  experimental  studies  by  Roever  and  Vining 
to  a  two-layered,  solid  half-space  covered  by  a  homogeneous  liquid 
cornterpart  [104-108]  or  by  a  liquid  layer  with  an  adjacent  gaseous  half¬ 
space  [109-110],  To  compose  the  theoretical  seismograms,  their  calcula¬ 
tions  took  into  account  up  to  two  subsequent  interactions  of  the  primarily 
radiated  P-wave  front  with  the  crucial  liquid/solid  interface.  Using  four, 
different,  small-scale  models  and  a  spark  gap  they  recorded  numerous 
pressure  responses  at  the  liquid/solid  interface.  In  almost  all  cases  the 
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measured  signals  again  showed  very  good  agreement  with  the  synthesized 
pressure/time  curves.  Usually  the  Scholte-wave  pulse  turns  out  to  be  the 
most  prominent  portion  of  the  signal  after  the  above-mentioned  direct 
arrival  in  the  liquid. 

The  extension  by  McLeroy  and  other  underwater  acousticians  of  these 
successful  laboratory  tests  to  the  more  complicated  and  contradictory  field 
experiments  in  shallow  water  will  be  described  in  Sect.  2.4. 
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2  EXPERIMENTAL  ASPECTS  OF  SEISMIC  SENSING  AT  THE  OCEAN  BOTTOM 

2.1  Geophones  as  Standard  Transducers  for  the  Detection 
of  Seismic  Signals 

The  classical  mass-spring  geophone  is  a  relatively  simple  electromechanical 
transducer.  Its  mechanical  system  consists  of  a  mass  that  is  suspended  by 
a  system  of  springs  within  a  frame  or  a  container.  For  different  technical 
reasons  the  motion  of  the  mass  has  to  be  restricted  by  an  axial  guidance, 
which  is  why  the  geophone  acts  always  as  an  unidirectional  sensor.  The 
amplitudes  of  this  axial  motion  have  also,  of  course,  to  be  limited  by 
damping  stops.  Usually  the  superimposed  electromagnetic  system  consists  of 
permanent  magnets  that  are  attached  to  the  housing  and  a  coil  of  wire  that 
is  wound  around  the  moving  mass.  The  magnets  are  fitted  to  the  frame  in 
such  a  way  that  they  create  optimal  field  conditions  around  the  moving 
coil.  As  the  moving  components  of  the  geophone  are  coupled  to  the 
container  both  by  the  springs  and  by  a  damping  dashpot,  its  response  curve 
is  characterized  by  a  flat  branch  for  all  higher  frequencies  and  by  a 
mechanical  resonance  for  a  relatively  low  frequency.  This  resonance  is 
usually  suppressed  partially  or  totally  by  a  suitable  electrical  load  but 
nevertheless  it  determines  the  low-frequency  cut-off  of  the  instrument. 

The  response  curve  of  the  sensor  decreases  for  all  frequencies  falling 
below  this  cut-off,  usually  at  vhe  relatively  high  rate  of  18  to  20 
dB/octave. 

Roughly  speaking  this  resonance  frequency  or  "natural  frequency"  is 
between  0.05  and  1  Hz  if  a  stationary  geophone  is  used  for  earthquake 
monitoring,  and  between  I  and  30  Hz  if  the  transducer  is  employed  for 
seismic  prospecting  or  surveillance.  By  putting  such  a  geophone  or 
seismometer  on  or  in  the  earth  we  try  to  create  a  rigid  connection  between 
the  housing  and  the  moving  environment.  Owing  to  its  inertia,  the 
suspended  mass  of  the  sensor  represents  the  "fixed  point  in  space"  relative 
to  which  we  want  to  determine  the  earth's  motion.  Due  to  the  above- 
mentioned  coupling  of  mass  and  container,  such  a  fixed  point  can  be 
established  only  approximately.  The  electrical  signal  at  the  output  of  the 
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coil  is  proportional  to  the  rate  of  relative  motion  of  magnet  and  coil  - 
thus  the  classical  geophone  is  primarily  a  velocity  detector.  Sometimes 
this  output  signal  is  immediately  differentiated  or  integrated  by 
incorporated  electronic  circuits  to  tune  the  sensor  to  acceleration,  or 
displacement  measurements.  The  interested  reader  may  find  the  pertinent 
mathematics  in  such  books  on  applied  geophysics  as  [111,112],  special 
papers  [113],  or  technical  publications  of  the  seismic  industry  [114] . 

As  seen,  the  geophone  acts  as  a  unidirectional  motion  detector  sensing  only 
the  velocity  component  parallel  to  its  axis.  The  resulting  sensitivity 
pattern  is  therefore  proportional  to  the  cosine  of  the  angle  measured  from 
the  axis.  Thus  it  has  a  dumbbell-like  shape  with  no  gain  in  any  direction 
but  with  increasing  attenuation  towards  the  off-axis  directions.  For  an 
omnidirectional  background  of  ambient  noise  and  an  on-axis  signal  field, 
the  signal-to-noise  ratio  has  been  calculated  to  be  approximately  5  dB 
higher  than  for  such  omnidirectional  sensors  as  most  of  the  commonly-used 
hydrophones  [115].  In  practice,  however,  this  signal-to-noise  ratio  will 
always  be  smaller,  because  the  sensitivity  pattern  of  real  geophones  never 
disappears  completely  in  the  off-axis  directions.  This  imperfection  is 
simply  a  consequence  of  a  slight  cross-axis  coupling  due  to  the  suspension 
and  axial  guidance  of  the  moving  mass  by  a  mechanical  spring  system. 

This  two-fold  function  of  the  springs  is  responsible  for  other  typical 
features  of  the  classical  geophone:  the  greater  delicacy  required  in  its 
handling  and  the  care  that  has  to  be  taken  in  adjusting  its  horizontal 
sensors.  Even  with  relatively  large  masses,  the  vertical  geophones  are 
quite  rugged  and  operational  even  for  alignments  between  ±5°  and  ±20°  from 
the  vertical.  On  the  other  hand  their  horizontal  counterparts  have  to  be 
handled  more  carefully,  be  positioned  correctly,  and  adjusted  within  a  few 
degrees. 

The  tilt  sensitivity  of  all  traditional  sensor  models  decreases  if  the 
resonant  frequency  is  increased  by  reducing  the  moving  mass  and/or 
stiffening  the  springs.  As  this  improves  the  general  ruggedness,  we  can 
easily  understand  the  success  of  short-period  phones  with  high  natural 
frequencies  (above  about  5  Hz)  in  the  field  of  seismir  prospecting.  It  is 
equally  obvious  why  the  use  of  high-performance  1  Hz  seismometers  has  been 
restricted  more  or  less  to  scientific  purposes. 
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In  1977  Teledyne-Geotech  marketed  a  1 -Hz  geophone  (model  S-500)  that  can  be 
transported  carefully  without  mass-locking  and  can  be  operated  in  any 
wanted  direction.  In  tne  past  such  favourable  characteristics  have  been 
offered  only  by  the  much  less  sensitive  piezoelectric  acceleration  sensors. 
With  these  new,  high-performance,  1 -Hz  geophones  it  is  now  possible  to 
reach  intrinsic  sensitivities  of  the  order  of  103  V/ m/s  by  using  closer 
winding  and  integrated  electronics. 

By  using  these  sensors  together  with  a  few  low-noise  amplifier  stages,  the 
seismic-noise  background  in  the  infrasonic  frequency  band  can  easily  be 
recorded  on  magnetic  tape  even  at  very  "quiet"  sites  on  land  or  on  the  deep 
ocean  floor  (see  Sect.  2.3)  where  the  vertical  deflection  (p-p)  lies 
between  approximately  1  my  at  1  Hz  and  0.01  my  at  30  Hz. 

Other  very  important  experimental  aspects  of  ocean-bottom  seismometers  are 
their  suitable  mechanical  design  (no  infrasonic  system  resonances),  their 
proper  coupling  to  the  sea  bed,  and  their  protection  against  spurious 
motion  effects  from  the  overlying  water  mass  [112,115-119] .  The  ideal 
condition  would  be  approached  by  using  very  small  embedded  sensors  that 
have  a  density  equal  to  the  average  of  the  surrounding  bottom.  Thus  the 
impedance  notching  should  be  optimized  by  a  proper  choice  of  weight  and 
shape  (especially  footpad  size)  of  the  receiver-station,  and  it  should  be 
as  well  sheltered  as  possible  from  currents  and  pressure  variations  in  the 
water  column  (e.g.  by  covering  with  acoustically-insulating  material). 

By  using  several  of  these  sensor  stations  on  or  just  below  the  sea  floor 
one  can  install  different  types  of  arrays  or  distributed  systems  and  thus 
profit  from  the  more  sophisticated  signal-processing  methods  (such  as 
beamforming)  developed  for  those  more  elaborate  receiver  configurations 
(see  e.g.  [115,116,). 


2.2  Some  Essential  Features  of  Sea-Floor-  Acoustics 

For  many  years  there  has  been  intense  research  on  sound  propagation  in  the 
ocean  bottom  and  the  literature  is  rich  with  reports  of  the  great  many 
theoretical  studies,  modelling,  and  field  experiments  performed.  We 
recommend  newcomers  to  the  field  to  the  summary  papers  by  Nafe  et  al 
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[120,121]  or  Hampton  [122],  and  the  more  modern  reviews  by  Hamilton 
[123,124]  or  Hampton  [125]  as  starting  points.  However,  because  of  the 
large  variety  of  marine  sediments,  and  the  unavoidable  deterioration  of 
model  experiments  (such  as  scaling  problems)  and  of  sample  tests  in  the 
laboratory  (using  piston-cores  or  grab-samples),  together  with  the 
technical  difficulties  of  in-situ  measurements,  there  still  remains  a 
substantial  lack  of  knowledge  about  the  subject,  particularly  within  the 
infrasonic  frequency  range  that  is  of  interest  in  the  present  context. 

We  therefore  do  not  want  to  encumber  the  reader  with  a  detailed  discussion 
of  the  numerous  theoretical  models  (such  as  [124,126,127 ,128])  and 
experimental  techniques  or  with  a  presentation  of  extensive  data  sets  (such 
as  [129,130]).  Thus  the  following  comments  are  intended  merely  to  give  a 
general  overview  of  the  main  features  and  some  specific  hints.  Speaking 
generally,  the  sea  floor  consists  of  a  series  of  "soft"  sediment  strata 
overlying  a  "hard"  rock  sub-bottom.  Under  real-world  conditions  this  scale 
of  hardness  is  of  course  very  smooth  and  the  transition  zones  may  be 
described  as  "relatively-hard"  sediments  or  "quite-soft"  rocks.  In  order 
of  increasing  depth,  the  lithosphere  below  the  ocean  may  be  classified 
roughly  as  follows: 

A.  water-saturated,  unconsolidated  sediment  layer(s)  (e.g.  mud,  silt, 
sand,  clay). 

B.  water-saturated,  consolidated  sediment  layer(s)  (e.g.  silt,  sand, 
clay). 

C.  sedimentary  rock  layer(s)  (e.g.  sand  or  limestone). 

D.  crystalline  rock  basement  (e.g.  basalt  or  granite). 

The  thicknesses  of  the  layers  in  A  and  B  may  vary  from  some  centimetres  to 
some  hundred  metres  (see  e.g.  [131])  while  the  primitive  rock  basement  D 
can  be  regarded  as  an  almost  homogeneous  half-space.  The  sedimentary-rock 
strata  of  group  C  lie  between  these  extremes.  The  consistency  of  these 
layer  materials  ranges  from  the  quasi-liquid  character  of  suspensions 
("very  soft")  to  the  extreme  hardness  of  crystals  ("extremely  hard").  For 
the  layers  in  A,  B  and  C  this  increasing  consolidation  is  mainly  caused  by 
the  weight  of  the  overburden  but  is  also  partly  due  to  chemical  processes. 
With  increasing  pressure,  the  packing  density  of  their  mineral  components 
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(with  many  different  grain  sizes  and  mixture  ratios)  becomes  greater  and 

consequently  the  contact  between  the  particles  becomes  more  rigid.  Thus 

the  resulting  or  mean  mass-density  p  of  the  layers  grows  from  about 

1.0  g/cm3  (water)  up  to  about  3.0  g/cm3  (crystalline  rocks),  the  porosity 

(percentage  of  porous  space)  dropping  from  almost  100%  (at  the  water/bottom 

interface)  to  0%  (in  sedimentary  rocks),  and  the  loose  contact  of  the 

particles  changes  into  a  quite  regular  solid  frame  or  internal  skeleton. 

From  this  it  is  evident  that  all  elastic  modulii,  and  thus  the  sound  speeds 

c  ,  c  ,  must  have  a  positive  gradient  with  growing  depth,  while  the  damping 
P  5 

losses  (mainly  due  to  internal  friction)  must  show  the  opposite  trend. 

These  very  general  considerations  of  the  acoustic  features  indicate  that 
the  ocean  floor  always  represents  a  "very  hard"  solid  if  we  penetrate 
sufficiently  deeply  and  may  be  composed  of  "quite  soft"  or  even  liquid-like 
media  near  the  water/bottom  interface.  If  we  disregard  the  influence  of 
pronounced  interfaces  and  inhomogeneities,  these  general  features  explain 
clearly  the  steady  upward  refraction  of  seismic  waves  and  the  often 
unexpectedly  good  conditions  for  long-range  transmission  via  deep  refracted 
paths.  Our  knowledge  of  the  acoustic  properties  of  the  deep  layers 
involved  originates  mainly  from  the  well-known  techniques  of  reflection  and 
refraction  shooting  [132  to  136]. 

As  here  we  are  mainly  interested  in  the  propagation  of  infrasound  at  or 
just  below  the  water/bottom  interface  we  will  concentrate  our  attention  on 
the  acoustics  of  the  upper  sea-floor  strata.  The  exploration  of  these  top 
layers  is  carried  out  by  means  of  the  already-mentioned  remote  techniques 
(see  also  [137-140]),  by  in-situ  implantation  of  probes  ([141-145]),  by 
laboratory  tests  of  grab  or  core-samples  (see  e.g.  [138, 140]),  or  by 
combined  methods  [146,147], 

At  this  point  we  would  like  to  offer  a  rough  idea  of  the  parameter  values 

that  govern  the  propagation  of  P-  and  S-waves  in  the  sea-floor.  To 

simplify  matters,  we  assume  a  normalized  sound  speed,  c  =  1500  m/s  at  the 

bottom  of  the  water  column,  and  relate  the  compressional  and  shear  velocity 

of  the  bottom  layers  to  it  by  two  constants  n  and  n  : 

P  5 

Cn  =  nn  Cw  i  =  n«  Cw  *  24> 
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For  the  same  reason,  we  generalize  the  well-known  high-frequency  behaviour 
of  many  materials  by  assuming  that  the  attenuation  coefficients  a  and 
(dB/m)  always  increase  linearly  with  frequency,  f  (kHz),  according  to: 

“p  *  V  ;  “s  =  Ksf-  <E<|-  26) 

Above  1  kHz  this  simple  relation  is  generally  valid  for  most  sedimentary 
bottom  strata,  as  the  numerous  data-points  of  Fig.  9  [127]  demonstrate.  In 
this  way  we  have  two  frequency-independent  attenuation  factors  K  and  K 

r  15 

(dB/km/Hz)  by  which  to  characterize  the  different  media.  Table  3  summarizes 
the  ranges  of  n  ,  Kp,  n$,  Kg  and  p  within  the  general  A-D  classification 
previously  described. 


TABLE  3  DENSITIES  p,  VELOCITY  RATIOS  n,  AND  ATTENUATION  FACTORS  K, 

FOR  TYPICAL  BOTTOM  LAYERS 

Compared  with  the  quite  steady  and  moderate  increase  of  density  with  depth 
[148),  the  values  of  the  other  parameters  increase  much  less  regularly  and 
cover  greater  ranges,  sometimes  even  of  orders  of  magnitude.  Because  of 
the  great  diversity  of  bottom  materials  the  limits  of  the  intervals 
indicated  in  Table  3  can  vary  widely,  so  the  figures  given  should  therefore 
be  regarded  as  indicative. 

The  mean  value  of  the  compressional  velocity  gradient  has  been  determined 
to  about  1  (s~^ )  in  the  sedimentary  layers  A-B-(C)  (to  a  depth  of  about  1  to 
2  km)  and  becomes  much  smaller  in  the  rock  basement  (C)-D  (see  e.g. 

[149> lt,0l) .  Thus  the  depth  profile  of  cp  is  much  steeper  than  that 
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of  the  overlying  water  column  and  explains  the  often-observed  strong  upward 
refraction  of  P-waves,  which  may  sometimes  even  create  apparent  low,  or 
even  negative,  bottom-reflection  losses  (see  e.g.  [151-153]).  The  low 
starting  value  of  0.97  for  np  demonstrates  a  well-known  suspension  feature: 
that  the  "disorder"  of  the  dispersed,  "fast",  solid  particles  in  a  "slow" 
liquid  may  even  decrease  the  original  sound  speed  of  the  dispersing  agent. 

The  dashed  line  in  Fig.  9  shows  a  me«n  value  of  Kp«  0.25  dB/km/Hz  for  the 
attenuation  of  P  waves  in  sediments.  With  growing  depth  the  attenuation 
factor  Kp  usually  decreases  quite  rapidly,  as  shown  in  Fig.  10,  and  then 
reaches  the  characteristic  value  for  rocks,  Kp  «  0.034  dB/km/Hz,  at  a  depth 
of  approximately  1.5  to  2  km  [154].  This  behaviour  is  a  consequence  of  the 
fact  that  the  total  loss  mechanism  is  more  or  less  reduced  to  the  basic 
material  losses  (internal  friction  and  scattering  effects)  because  those 
due  to  the  friction  of  the  grains  in  a  sort  of  "low-frequency  skeleton" 
("external"  friction)  and  viscosity  are  eliminated  stepwise  by  compaction. 
Table  3  gives  only  one  Kp  interval  for  all  rock  types,  because  it  may 
sometimes  happen  that  intermediate  sedimentary  rock  layers  are  even  less 
lossy  than  their  crystalline  basement.  More  modern  studies  [128]  call  in 
question  the  validity  of  Eq.  25  for  the  infrasonic  frequency  regime  and 
thus  assume  much  smaller  Kp  values  in  the  "soft"  top-layers  due  to  the 
above  mentioned  "low-frequency  skeleton"  which  is  thought  to  be  a  macro¬ 
scopic  counterpart  to  the  microseismic  lattice  structure  of  crystals. 

The  ns~values  in  Table  3  indicate  that  the  shear  velocity  in  unconsolidated 
sediments  may  be  very  low  (down  to  about  5%  of  the  sound  speed  in  the 
overlying  water),  but  even  the  extremely  small  amount  of  data  concerning 
the  depth  dependence  suggests  the  already  expected  high  velocity  gradient 
of  about  10(s’b  for  the  first  10  m  and  of  about  3(s"^)  for  the  following 
100  m  [155,156].  The  often  neglected  existence  of  S-waves  in  those  very 
"soft"  top-layer  sediments  is  usually  explained  with  help  of  a  sort  of 
macroscopic  shear  rigidity  assigned  to  the  "low-frequency  skeleton". 

However,  the  pertinent  attenuation  is  very  high,  as  Table  3  and  the  few 
published  data  sets  [157]  show,  and  a  mean  value  given  by  Ks«ll  dB/km/Hz 
may  be  quite  reasonable  ir  this  case.  With  growing  depth,  the  attenuation 
factor  Kg  is  also  expected  to  decrease  extremely  fast,  because  we  know 
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FIG.  10  DEPTH  DEPENDENCE  OF  THE  ATTENUATION  FACTOR  Kp  FOR  COMPRESSIONAL 
waves  in  sediments  (data  compiled  from  literature) 
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(e.g.  from  cores  taken  during  the  Deep  Sea  Drilling  Project)  that  in  the 
underlying  rock  materials  Kg  adopts  values  of  the  same  order  as  Kp. 

As  the  top  layers  of  the  sea  floor  usually  belong  to  the  stratifications  A, 

B  and  C,  the  velocities  of  the  infrasonic  interface  waves  ("pure"  or 
modified  Scholte  waves;  see  Table  1)  may  range  from  about  csc^  «  60  m/s  in 
case  of  relatively  thick  (at  least  more  than  one  wavelength)  and  "very- 
soft"  sediment  covers  (c$«cw),  to  that  of  the  sound  speed  in  water 
csch  ~  1500  m^s  in  case  "quite  hard"  sedimentary  rocks  (cg  >  cw). 

Because  they  result  from  a  combined  action  of  compression  and  shear,  these 
interface  waves  are  governed  by  the  usually  quite  high  shear  attenuation. 
Thus  we  cannot  expect  long-range  or  even  medium-range  propagation  of  this 
wave  type  in  areas  where  the  basement  is  covered  by  very  thick  and 
"extremely  soft",  unconsolidated  sediment  strata  [147]. 

In  case  of  a  "hard"  rocky  sea  floor  that  may  even  have  a  thin  and  very 
lossy  sediment  cover  ("coating"),  the  propagation  conditions  are  certainly 
much  better  and  may  allow  the  successful  use  of  those  guided  waves  as  a 
tool  for  the  seismic  sensing  of  waterborne  infrasound,  (see  Sect.  3.2). 

Unfortunately  these  more  promising  aspects  of  a  rock-like  bottom  are  merely 
based  on  the  material  parameters  and  do  not  account  for  the  macroscopic 
structure  and  shape.  Under  real-world  conditions  such  a  hard  sea  floor  is 
not  only  restricted  to  relatively  small  areas  (e.g.  banks  and  reefs)  but  is 
mostly  made  up  of  inclined  and  sometimes  irregular  layers  and  therefore 
presents  a  surface  full  of  features.  These  latter  features  are  particularly 
obstructive  to  the  excitation  and  propagation  of  all  guided  seismic  waves. 
Thus  the  optimal  conditions  for  a  successful  excitation  and  propagation  of 
interface  waves  may  be  found  in  those  wide-spread  areas  where  the  marginal 
zone  of  the  sea  floor  is  composed  of  a  quite  regular  sediment  layering 
providing  a  smooth  surface  and  a  rapid  consolidation  with  growing  depth. 
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2.3  Acoustic  and  Seismic  Ambient  Noise  Data 
from  Shallow-Water  Areas 

* 

In  sensing  acoustic  or  seismic  signals  we  must  always  operate  against  a 
*  certain  noise  background  that  is  of  many  origins  (sources),  qualities 

(spectra),  and  quantities  (levels).  As  with  sea-floor  acoustics,  the  study 
of  ambient-noise  characteristics  on  land  and  at  sea  is  a  large  field  and  we 
must  restrict  ourselves  to  some  highlights  and  to  referring  the  reader  to 
the  pertinent  literature. 

Since  World  War  II  underwater  acousticians  have  devoted  much  effort  to  the 
registration  and  interpretation  of  waterborne  or  acoustic  noise,  and 
seismologists  have  acted  likewise  with  regard  to  earth  movements.  As 
summary  papers  on  the  first  subject  we  recommend  the  book  by  Ross  11581, 
the  article  by  Arase  11591  and  the  classical  paper  by  Wenz  11601  .  A  good 
insight  into  the  nature  of  seismic  noise  spectra  is  given  by  the  publica¬ 
tions  of  Frantti  et  al  1161,1621  and  Brune  et  al  11631.  Both  the  spatial 
division  of  these  research  domains  and  the  interest  in  different  frequency 
r  bands  caused  quite  poor  overlapping  of  the  data  within  the  range  between 

1  Hz  and  100  Hz.  With  a  growing  interest  of  the  sonar  community  in  the 
low-frequency  regime  and  the  development  of  more  versatile  seismic  sensor 
stations  (OBS's)  this  gap  started  to  be  filled  in  the  early  1960s.  While 
the  underwater  acousticians  extended  their  noise  measurements  into  coastal 
waters  [164-166]  the  seismologists  have  been  mainly  concerned  with  data 
collections  in  the  deep  ocean  floor  1167-1711 . 

Consequently  our  knowledge  of  infrasonic  seismic  noise  in  shallow  water  is 
still  quite  poor  and  based  on  a  few  field  trials.  Another  shortcoming  of 
the  existing  seismic  data  is  given  by  the  fact  that  most  of  them  have  been 
recorded  merely  with  vertical  geophones  (see  Sect.  2.1). 

We  will  now  reproduce  a  few  representative  curves  and  discuss  their  most 
prominent  features.  Figure  11,  taken  from  [166]  ,  recapitulates  the  well- 
known  pressure  spectrum  levels  of  the  infrasonic  and  low-frequency  noise  in 
shallow  water  and  indicates  the  main  noise  sources. 

With  increasing  frequency  the  acoustic  noise  level  generally  drops  at  a 
t  rate  of  -8  to  -10  dB/oct.  Below  10  Hz  the  wind-generated  pressure 
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FIG.  11  SOUND-PRESSURE  SPECTRUM  LEVELS  FOR  DIFFERENT  NOISE  SOURCES 
IN  THE  SEA 
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fluctuations  (surface  waves  and  turbulences)  form  the  dominant  source 
mechanism,  while  above  10  Hz  the  traffic  noise  becomes  increasingly 
important.  As  a  seismic  counterpart,  we  have  plotted  in  Fig.  12  two  pairs 
of  curves  and  S2  based  on  measurements  by  McLeroy  [172]  at  two  different 
shallow-water  locations.  The  smooth,  dashed  curves  give  the  averaged 
velocity  spectrum  levels  of  the  vertical  sensor  while  the  undulating,  solid 
curves  represent  (two  of  the  few  examples  for)  the  averaged  horizontal 
levels.  The  magnitude  of  the  horizontal  components  exceeds  that  of  the 
vertical  ones  by  approximately  10  dB,  which  is  also  typical  for  propagation 
studies  (see  Sect.  2.4)  and  may  be  a  consequence  of  trapped  shear  waves  in 
sediments  [173] .  However,  the  "snapshot"  character  of  these  curves  does 
not  permit  far-reaching  conclusions.  With  increasing  frequency  these 
seismic  noise  levels  indicate  a  fall -off  that  amounts  to  only  about 
-6  dB/oct  and  thus  seems  to  be  less  than  in  the  acoustic  case.  Figure  13 
presents  another  simplified  plot  [172]  that  compares  the  average  S  of  the 
dashed  curves  in  Fig.  12  with  different  deep  sea  results  D^,  and  with 
average,  T,  and  minimum  values,  LQ,  on  land.  Curve  is  based  on  data 
from  [171] ,  Do  from  [169] ,  Dj  from  [168],  and  from  [167],  and  T 
and  Lq  from  [163] .  As  the  most  important  noise  sources  that  concern  the 
deep  ocean  bottom  are  formed  by  microseisms  (mainly  due  to  moving  water 
masses),  by  biological  activity,  and  by  distant  shipping,  the  levels  there 
are  not  very  different  from  those  on  land,  merely  that  the  fall -off  with 
increasing  frequency  is  much  steeper  (approximately  -10  to  -20  dB/oct)  than 
it  is  for  land  sites  (about  -8  dB/oct).  This  latter  discrepancy  did  not 
occur  in  connection  with  the  Black  Sea  data  published  by  Rykunov  et  al 
[ 170]  for  the  band  from  2  to  15  Hz. 

Bradner  ( 1741  has  shown  that  the  feed-back  from  these  microseisms  creates 
the  often-observed  increase  in  the  acoustic  noise-level  towards  the 
water/bottom  interface  (positive  gradient  of  the  depth  profile)  due  to  the 
microseismic  Scholte  waves. 

The  spectrum  levels  of  seismic  noise  in  coastal  waters  are  expected  to  be 
higher  than  those  measured  at  the  deep  ocean  floor.  The  shallow  bottom 
regions  are  much  closer  to  the  moving  sea  surface,  to  the  routes  of  heavy 
ship-traffic,  to  offshore  platforms,  and  to  the  shore  line  with  its  natural 
unrest  (e.g.  breaking  waves)  and  its  man-made  noise  sources  (industrial  and 
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FREQUENCY ,  Hz 


FIG.  12  AVERAGED  VERTICAL  AND  HORIZONTAL  VELOCITY  SPECTRUM  LEVELS  ON 
THE  SEA  FLOOR  IN  SHALLOW  WATER  [172] 


FREQUENCY. Hz 


FIG.  13  AVERAGED  VERTICAL  VELOCITY  SPECTRUM  LEVELS  FOR  DIFFERENT  SITES 
ON  LAND  AND  ON  THE  SEA  FLOOR  IN  SHALLOW  AND  DEEP  WATER  [172] 
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cultural  activities).  This  trend  was  also  confirmed  by  two  more  shallow 
sensor  installations  (at  500  and  300  m  water  depth)  during  the  above- 
mentioned  Black  Sea  experiment  (the  levels  were  about  20  dB  above  the 
deep- bottom  results).  More  actual  deep-ocean  data  published  by  Asada  et 
al  [175]  confirm  the  steep  decay  of  the  spectrum  level  in  the  range  from  1 
to  10  Hz  (about  -25  dB/oct)  but  their  magnitudes  are  almost  30  dB  higher 
and  thus  exceed  even  the  shallow-water  records  presented  above.  In 
addition,  their  curve  passes  a  pronounced  minimum  at  about  15  Hz  and  then 
again  shows  a  modest  increase  towards  higher  frequencies.  Some  preliminary 
results  from  noise  measurements  with  a  prototype  of  SACLANTCEN' s  ocean- 
bottom  seismometer  in  Italian  coastal  waters  (at  about  20  m  depth  off 
Viareggio)  seem  to  confirm  the  existence  of  such  a  minimum  in  the  seismic 
noise  level  between  10  Hz  and  20  Hz. 

The  acoustic  counterpart  of  this  noise  minimum  has  been  found  both  at 
higher  (at  about  25  Hz  in  [164])  and  at  lower  frequencies  (at  about  6  Hz  in 
[165]). 

Besides  the  favourable  aspects  of  seismic  paths  for  infrasound  propagation 
(see  Sect.  3.2)  this  "low-noise  window"  (just  in  the  middle  of  the 
frequency  band  of  interest)  may  play  a  decisive  role  in  the  successful 
application  of  seismic  sensing  to  detection  purposes. 


2.4  Field  Experiments  Concerning  Seismic  Detection 
of  Waterborne  Infrasound 


The  results  of  some  of  the  classical  shallow-water  studies  [70,74,85]  and 
developments  in  the  detection  of  microseisms  (see  Sect.  2.3)  and  nuclear 
blasts  [1761  stimulated  a  few  underwater  acousticians  to  make  various  small 
tests  and  to  perform  some  quite  extensive  experiments  under  real-world 
conditions  to  explore  low-frequency  transmission  in  shallow  water. 
Unfortunately,  in  its  early  days  this  field  work  had  quite  different 
objectives  from  the  present  subject,  often  put  no  efforts  on  the  collection 
of  environmental  data,  and  partially  suffered  from  inadequate  equipment. 

The  latter  technical  shortcoming  refers  especially  to  the  horizontal 
geophones,  which  frequently  had  to  be  excluded  because  of  the  need  for 
delicate  handling  (see  Sect.  2.1). 
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During  the  years  1943/44  Worzel  and  Ewing  [17?]  installed  a  sensor  station 
with  two  hydrophones  and  one  vertical  geophone  on  the  ocean  bottom  in 
different  shallow-water  areas.  At  water  depths  between  15  m  and  180  m  they 
sensed  the  pressure  and  displacement  signals  from  explosive  charges  dropped 
at  distances  of  between  0.5  and  50  km.  Their  main  interest  was  devoted  to 
the  so-called  "ground-wave"  (85] ,  which  is  composed  of  the  shallow- 
refracted  arrivals  having  dominant  frequencies  in  the  10  to  30  Hz  band. 

About  1958  Worley  and  Walker  1178]  performed  a  similar  study  with  a  more 
versatile  receiver  package  but  unfortunately  their  horizontal  geophone 
failed.  Nevertheless  they  were  able  to  show  that  the  signals  from  the 
hydrophone  and  the  vertical  geophone  have  the  same  range-dependent  fall -off 
in  the  15  to  70  Hz  frequency  regime. 

Approximately  at  the  same  time,  Blaik  and  Clay  [179]  tried  to  sense 
waterborne  sound  on  land.  They  towed  powerful  low-frequency  transducers  in 
shallow  water  and  installed  vertical  geophones  in  a  well  on  shore  at  depths 
between  about  100  and  600  m.  Profiting  from  the  refracted  field  at  shallow 
and  medium  depth  they  succeeded  in  detecting  CW  signals  (10,  24,  88  and 
148  Hz)  over  ranges  up  to  about  9  km. 

From  1960  to  1963  Shorthouse  l 1801  developed  and  successfully  applied  a 

seismic  system  to  observe  P-head  waves  from  velocity  discontinuities  in  the 

oceanic  crust  and  layering  in  the  sediment  overburden.  Using  a  geophone 

and  several  hydrophones  on  the  deep  sea  floor  he  recorded  the  signals  from 

6  lb  (2.3  kg)  charges  exploded  on  the  ocean  bottom  at  different  distances. 

As  a  by-product,  long,  late-arri ving  wave  trains  were  observed  up  to  ranges 

of  2.5  km;  these  showed  a  substantial  dispersion  and  quite  strong  amplitude 

-1  4 

decay  proportional  to  r  '  .  Analyzing  the  pertinent  data,  Davies  [181] 
clearly  identified  these  "tails"  of  the  records  as  Scholte  waves  with  phase 
velocities  between  130  and  60  m/s  in  the  frequency  interval  from  2  to 
10  Hz. 

Since  1958  McLeroy  and  co-workers  have  carried  out  different  theoretical 
studies,  model  experiments  (see  Sect.  1.5),  and  quite  extensive  field  tests 
on  use  of  seismic  paths  for  detection  and  communication  in  coastal  waters. 
Using  a  shallow,  buried,  triaxial  0BS  station  and  a  nearby  hydrophone  they 
recorded  the  signals  from  high-power  broadband  CW-sources  and  explosives. 
Despite  the  fact  that  their  model  studies  successfully  synthesized  the 
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different  contributions  to  the  sound  field  at  an  interface,  they  were  often 
unable  to  ascribe  the  complex  real-world  data  to  these  well-defined  wave 
types.  In  general,  they  were  therefore  merely  able  to  distinguish  between 
the  ensembles  of  the  "ground-wave"  and  the  "water  wave".  However,  they 
noticed  that  at  all  frequencies  there  was  a  much  higher  signal-to-noise 
ratio  on  the  horizontal  geophones  (up  to  20  dB)  than  on  the  two  other 
sensors.  During  one  of  these  field  experiments  they  observed  a  distinct 
interface  wave  that  they  interpreted  as  a  Scholte  wave  despite  some 
deviations  from  the  expected  characteristics  [182]  .  The  pertinent  wavelets 
occurred  at  frequencies  between  6  and  7  Hz  and  at  short  distances  (up  to 
600  m)  their  propagation  loss  was  slightly  smaller  (about  3  dB)  than  that 
of  the  sound  energy  being  trapped  in  the  water  column.  This  latter  feature 
is  in  contradiction  to  some  other  sea-floor  studies  that  used  the  interface 
wave  as  a  probe  and  has  disclosed  very  high  losses  even  for  the  infrasonic 
frequency  regime  (see  Sect.  2.2). 

In  1S68  Urick  1183,184]  tried  to  track  small  craft  (controlled  runs)  and 
large  freighters  (uncontrolled  runs)  with  a  similar  sensor  package  on  the 
sea  floor  (at  about  18  m  depth)  in  an  area  with  a  rapidly  sloping  bottom. 

For  all  frequencies  below  100  Hz  the  signal-to-noise  ratio  of  the  vertical 
geophone  was  less  than  that  of  the  hydrophone  and  showed  a  very  rapid  fall- 
off  with  range,  while  that  of  horizontal  sensors  was  always  greater  and 
decreased  more  slowly  with  range.  None  of  the  results  so  far  mentioned  has 
given  an  indication  of  strong  interface  waves  and  all  appear  bewildered  by 
the  unexpectedly  high  levels  of  the  cross-deflection  (up  to  50%  of  the 
radial  displacement  [1831).  The  latter  phenomenon  was  generally  explained 
by  a  relatively  high  energy  transfer  to  Love  waves  caused  by  inhomogeneities 
and  a  non-regular  layering  (see  Sect.  1.4). 

Approximately  at  the  same  time  0.  Hastrup  and  T.  Akal  of  SACLANTCEN  had 
already  made  a  systematic  test  in  the  Gulf  of  La  Spezia  to  investigate  the 
existence  of  an  explosion-generated  Scholte  wave  at  the  water-bottom 
interface.  Unfortunately  their  efforts  failed  because  of  the  shortcomings 
of  the  equipment,  which  consisted  of  quite  insensitive  vertical  sensors 
(prospecting  geophones)  in  conjunction  with  very  modest  display  and 
recording  facilities  (oscilloscopes  and  photographic  recorders). 
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Despite  our  devaluation  of  the  application  of  vertical  geophones,  it  must 
be  noted  that  Bucker  [185]  has  recently  demonstrated  that  such  bottom- 
mounted  sensors  are  well  suited  to  detect  the  infrasonic  lines  of  passing 
ships  in  shallow  water. 

In  1976  McLeroy  and  his  co-workers  [186]  again  performed  a  very  extensive 
real-world  experiment  in  the  Gulf  of  Mexico.  Using  a  receiving  array  of 
nine  triaxial  OBS  stations  and  six  bottom-mounted  hydrophones,  they 
conducted  numerous  acoustical  runs  in  shallow  water  with  towed,  powerful  CW 
sources  (operated  at  distances  up  to  about  50  km)  and  some  200  explosive 
charges  at  20  m  depth  and/or  1  m  above  the  sea-floor  (detonated  at 
distances  of  up  to  about  16  km).  This  time  the  environmental  background 
was  very  carefully  checked  by  taking  many  samples  of  bathymetry,  sound- 
speed  profile,  wind  speed,  wave  height,  currents,  bottom  layering,  and  sub¬ 
bottom  structure  (see  Sect.  2.  2). 

Until  the  end  of  1978  only  a  few  of  these  data  had  been  analyzed  but  the 
preliminary  results  already  indicated  a  muci.  better  signal-to-noise  ratio 
for  the  radial  and  the  vertical  geophone  than  for  the  transverse  sensor  and 
the  nearby  hydrophone.  This  contradiction  to  the  outcome  of  all  former 
field  tests  has  to  be  confirmed  by  the  remaining  data  and  investigated  in 
detail.  On  the  other  hand,  this  surprising  result  is  typical  of  the 
complexity  of  seismic  sensing  and  of  the  associated  experimental  diffi¬ 
culties. 

In  1976/77  a  simple  triaxial  OBS  prototype  was  developed  by  SACLANTCEN’s 
Environmental  Acoustics  Group  using  three  highly-sensitive  1 -Hz  geophones 
of  the  Geospace  HS-10/1B  type  and  one  nearby  hydrophone.  During  some 
shallow-water  trials  in  1977  and  1978  analogue  data  of  background  and  ship 
noise  as  well  as  signals  from  small  explosive  charges  were  transferred  via 
cable  or  radio  link  to  a  receiving  ship  [187].  Owing  to  many  problems  with 
the  very  delicate  sensors,  the  mechanical  design  of  the  station,  and  the 
electronics,  our  initial  experimental  results  have  not  been  satisfactory 
(see  Sect.  2.3) .* 


Since  wtiLUng  thcb  text,  thebe  ptcblemb  have  been  oveA.com  and  bucceAb^ul 
ieZd  woAk  mb  pet  foamed  in  1919  (bee  Appendix  A  de,tailb). 
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In  most  of  the  above-mentioned  field  trials  the  OBS  station  was  fitted  with 
a  hydrophone  to  provide  a  better  understanding  of  the  propagation  phenomena 
(clear  distinction  between  the  different  wave  types,  which  in  most 
instances  has  not  yet  been  achieved)  and/or  to  get  a  true  standard  of 
comparison  between  the  detection  methods  (direct  comparison  between  the 
acoustic  and  seismic  signal-to-noise  ratios). 

This  procedure  proved  successful  in  the  famous  Vela  Uniform  Program  (of 
the  Advanced  Research  Projects  Agency)  for  the  monitoring  of  large  under¬ 
ground  explosions  (see  also  [188]),  and  it  may  be  extremely  helpful  as  a 
means  of  eliminating  or  profiting  from  the  influence  of  converted  or 
trapped  shear  waves  (see  Sects.  L.  2  and  2.3). 
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3  MODELLING  OF  SEISMIC  WAVE  PROPAGATION  IN  THE  SEA-FLOOR 

3.1  Characteristics  of  the  Thomson-Haskell  Matrix  Method 

In  Ch.  I,  Sects.  1.1,  1.3  and  1.5  attempted  to  describe  the  very  complex 
seismic  field  as  a  whole  and  Sects.  1.2  and  1.4  concentrated  on  interface 
waves.  For  a  more  realistic  evaluation  of  the  role  of  different  wave  types 
and  more  reliable  predictions  for  field  experiments,  we  have  to  fit  all 
available  environmental  data  into  a  theoretical  model.  Such  a  model  should 
be  as  comprehensive  as  possible  and,  at  the  same  time,  well  suited  for 
computer  application. 

Section  1.3  suggested  that  the  separation  of  the  different  contributions  to 
the  sound  field  by  means  of  classical  mathematics  would  often  be  very 
elaborate,  sometimes  extremely  difficult,  and,  moreover,  quite  inexact  in 
most  cases.. 

If  we  treat  the  earth's  crust  as  a  horizontally-stratified  half  space 
consisting  of  homogeneous  liquid  and  solid  layers  with  constant  parameters 
(density  p,  sound  speeds  c  and  c  ,  attenuation  coefficients  a  and  a  ), 

r  **  r 

all  these  requirements  are  met  by  the  well-known  Thompson-Haskel 1  matrix 
method  [  7, 189,iu0) .  Here  we  will  give  merely  a  snort  outline  of  its  basic 
principles . 

Under  simple  geometrical  conditions  (plane  waves  incident  from  infinity,  or 
else  cylindrical  or  spherical  waves  radiated  by  a  horizontal  or  vertical 
line  source  or  by  a  point,  source)  the  solutions  to  the  wave  equations  in 
homogeneous,  source-free  layers  are  well  known  IS  to  8,  221.  The  classical 
way  to  solve  those  range  (x  or  r),  depth  (z),  ar.d  time  (t)-dependent 
propagation  problems  is  the  application  of  a  Fourier  tran'form  [SS-SS]  for  the 
range  and  a  Laplace  transform  (b<-56]  for  the  time  in  the  case  of  impulsive 
sources.  Thus  the  original  partial  differential  equations  are  converted  to 
algebraic  equations  with  respect  to  the  complex  wave  number  (k)  space  and 
the  complex  object  space  of  the  Laplace  transform.  There  the  resulting 
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ordinary  differentia]  equations  concerning  z  can  be  solved  easily  with 
regard  to  the  similarly  transformed  boundary  conditions  due  to  the 
separation  of  the  original  variables.  Most  of  the  mathematical  diffi¬ 
culties  usually  arise  from  the  reverse  transform  of  the  solution  back  to 
the  physical  space/time  domain,  and  therefore  classical  approaches  (71, 
computer-aided  numerical  methods  11911,  or  unusual  concepts  H27have  to  be 
used. 


As  the  harmonic  time -dependence  of  incident  waves  can  be  separated  and 
thereby  excluded  immediately,  Thomson  1189]  and  Haskell  11901  have 
demonstrated  by  means  of  monochromatic  plane  waves  that  the  resulting  sound 
field  within  the  m-th  layer  can  be  represented  by  a  four-vector 

in  the  k-space.  Because  it  arises  from  the  interference  of  obliquely- 
downward  and  upward-travelling  waves,  this  vector  may  be  related  to  the  z- 
dependent  amplitudes  of  the  basic  potentials  (displacement  or  velocity 
potential  of  P-  and  S-  waves)  or  to  such  intermediate  terms  as  the  z- 
dependent  amplitudes  of  the  two  particle  velocities  u ,  w  and  the  two 
stresses  xZ2  =  p,  t  =  T  according  to: 


/s 


\  Tm  I 


(Eq.  26) 


This  latter  form  of  the  layer  vector  is  obvious';  '  the  most  convenient 
one  because  it  already  contains  those  terms  that  have  i  be  continuous  at 

the  interfaces.  Being  composed  exclusively  from  the  solutions  to  a  system 

of  two  wave  equations,  the  vertical  "propagation"  through  the  pertinent 
homogeneous  layer  with  constant  parameters  is  also  well  defined.  Thomson 
and  Haskell  have  shown  that  this  correlation  implies  relatively  simple  and 
favourable  mathematical  consequences.  The  vector  at  the  bottom,  m,  of  the 
layer,  m,  can  be  calculated  from  its  counterpart  at  the  top,  m,  or 
viceversa  by  a  simple  matrix  multiplication  (linear  operation): 

i  =  jS  £-  ■  £  (Eq.  27) 

m  m  m  ’  m  n:  m 

The  4  x  4  layer  matrix  &  is  a  function  of  the  medium  parameters,  the 

layer  thickness  Hm,  the  complex  wave  number  j<,  and  the  resulting  phase 
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velocity  c,  and  thus  includes  the  effects  of  damping  and  dispersion.  The 
interested  reader  may  find  the  explicit  form  of  the  layer-matrix  elements 
a^v  in  the  repeatedly-cited  original  papers,  in  the  book  by  Bath  [7] ,  and  in  a 
report  by  Kutschale  [192] .  As  the  boundary  -  or  transition  -  conditions  at 
the  top  and  the  bottom  of  layer  m  require  the  fulfilment  of  the  relations: 


A  A  A  A 


(Eq.  28) 


and  Eq.  27  is  valid  for  all  three  adjacent  layers,  we  can  also  link  the 

A  A 

vector  ^  m+i  to  its  more  remote  counterpart  by  matrix  multiplication: 


£ 

m+1 


=  ,  A  Jl  , 

m+1  m  m-1 


=  A 


syst 
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syst 


£ 


m+1 


(Eq.  29) 


Thus  the  resulting  matrix  -^SySt  of  a  stack  of  layers  (each  of  them  with 
constant  parameters)  is  simply  the  product  of  all  section  matrices,  and  we 
can  finally  link  the  radiation  condition  (allowing  only  outgoing  waves)  of 
the  deepest  layer  (forming  the  sub-bottom  half-space)  to  the  pressure- 
release  surface  of  the  top  layer  (earth  or  sea  surface)  or  viceversa: 


= 

n 

&  syst 

(Eq. 

30a) 

a 

£  = 

-1  A 

A  1 

(Eq. 

30b) 

T 

syst  n 

Proceeding  in  this  way,  all  solution  elements  are  compatible  witn  all 
transition  conditions  and  we  thus  obtain  a  mathematical  description  (in  the 
complex  k-space)  of  the  complete  wave  field. 


The  above-mentioned  fact  that  each  layer-matrix  Am,  and  accordingly  the 
resulting  stack-matrix  is  making  allowance  for  dispersion  may  be 

easily  checked  by  considering  its  version  for  extremely-low  and  extremely- 
high  frequencies,  17,  190] .  Iri  the  low-frequency  limit  the  resulting  field 
is  identical  with  a  Rayleigh-type  wave  at  the  surface  of  a  layered  half¬ 
space  1 193 ]  and  in  the  high-frequency  limit  it  splits  into  the  typical 
interface  waves  at  the  various  interfaces. 
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Because  of  the  exclusion  of  range  dependence,  the  above-defined  field 
vectors  and  layer  matrices  can  be  used  in  connection  with  all  two- 
dimensional  wave-propagation  problems.  In  practical  operation  these 
mathematical  tools  have  to  be  applied  in  the  following  way:  the  surface- 

A  A 

vector  has  to  be  calculated  from  the  deepest  sub-bottom  vector  £■- 
according  to  Eq.  30b  by  "carrying"  the  field  vector  across  all  layers 
and  interfaces.  The  field  vector  then  has  to  be  "shifted"  down  again  to 
the  depth  at  which  the  seismic  signal  is  of  interest  by  using  a  reduced 
system  matrix  in  Eq.  30a.  Thus  it  is  obvious  that  the  receiver  depth  must 
coincide  with  the  coordinate  of  an  actually  existing  interface  or  of  an 
artificially  inserted  auxiliary  boundary  without  physical  meaning.  Similar 
arguments  are  valid  for  the  positioning  of  the  source  that  is  of  most 
practical  interest:  the  point  source  SQ  that  radiates  spherical  P-waves. 

In  the  range-independent  object  space  such  a  source  provokes  a  well-defined 
discontinuity  of  the  vertical  particle  velocity  w.  This  discontinuity  can 
be  incorporated  in  our  formalism  by  inserting  another  extra  interface  that 
divides  the  original  source-layer  s  into  two  sub-layers  s'  and  s", 
and  by  defining  a  pertinent  "source  vector"  P  to  be  added  to  the  field- 

A 

vector  £-„ 

l  -  £_  +  fQ  (Eq.  31) 

r  *  0 


during  the  composition  of  the  total  field  by  means  of  Eq.  30b.  In  the  case 
of  a  harmonic  point  source  this  source  vector  p  is  of  the  simple  form 

(861: 


Finally,  to  establish  the  encoded  range  dependence  of  the  resulting  wave 
field  (which  is  independent  of  the  bearing),  and  thus  to  describe  the 
seismic  signal  with  respect  to  a  remote  receiver  position,  we  have  to 
perform  an  inverse  Fourier  transform,  which  can  be  carried  out  using  modern 
computer  techniques  [191,192] . 
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3.2  Some  Preliminary  Results  from  the  Fast-Field-Program 

In  Sects.  1.4  and  2.2  we  have  already  mentioned  that  realistic  bottom 
conditions  usually  create  a  complex  eigenvalue  problem  (with  regards  to 
sound  propagation  in  the  shallow-water  channel)  that  can  not  be  solved  by 
the  classical  normal-mode  technique.  The  interesting  Scholte  wave  (zeroth 
mode)  always  corresponds  to  a  real  eigenvalue  but  neither  can  it  be  handled 
in  a  proper  way  because  of  basic  limitations  in  the  underlaying  theory. 
Fortunately  this  and  other  seismic  problems  can  be  tackled  with  the  help  of 
the  above-described  Thompson-Haskell  matrix  method.  Harkrider's  solution 
technique  of  the  wave  equation  (194]  and  Dorman's  matrix  notation  (1931 
simplified  the  application  of  this  many-membered  method  but  its  successful 
breakthrough  did  not  become  possible  until  the  introduction  of  the  Fast- 
Field-Program  technique  by  Marsh  and  Elam  (1951  and  by  Di  Napoli  (1961. 

At  this  point  we  would  like  to  demonstrate  the  capabilities  of  such  an  FFP, 
made  available  to  SACLANTCEN  by  Kutschale  [192]  from  Lamont-Doherty 
Geological  Observatory.  The  version  implemented  on  SACLANTCEN's  UNIVAC 
1106  computer  1197]  is  able  to  calculate  the  complete  monochromatic  wave 
field  in  a  stack  of  isovelocity  liquid  and  solid  layers  using  optional 
depths  of  source  and  receiver. 

To  work  out  some  of  the  essential  features  and  trends  of  the  interesting 

seismic  contribution  to  long-range  transmission  we  have  chosen  a  very 

simple  but  realistc  type  of  sea  bed  in  which  a  sand  layer  (p  =  2.0  g/cm3; 

c  =  1800  m/s;  c  =  700  m/s)  with  high  absorption  (a  =  0.75  dB/X; 

_ P  ^  P 

as  =  1.5  dB/X)  and  constant  thickness  =  L  rests  on  top  of  a  sedimentary 

rock  (p  =  2.2  g/cmJ;  Cp  =  3000  m/s;  c$  =  1400  m/s;  =  0.1  dB/X; 

=  0.05  dB/X;  =  ®).  We  assume  a  10  Hz  sound  source,  SQ,  in  the 

middle  or  at  the  bottom  of  a  water-column  (p  =  1.0  g/cm3,  c  =  1500  m/s; 

—  -6  ” 
a  -  1.8  x  10  dB/X)  of  mostly  =  W  =  60  m  depth,  with  the  receiver 

always  placed  on  the  ocean  floor. 

For  extreme  simplification,  all  wave  types  propagating  faster  than  the 
speed  of  sound  in  water  are  included  in  the  "acoustic  portion"  of  the  field 
while  all  slower  types  are  counted  as  the  "seismic  portion".  Such  a 
drastic  splitting  of  the  total  elastic  wave  field  can  be  justified  by  means 
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of  Fig.  7  where  the  phase  velocities  of  the  blurred  modes  in  the  water  are 
seen  to  be  well  separated  from  those  of  the  interface  wave.  The  present 
example  is  more  complex  of  course,  and  such  a  combination  obscures  the  fine 
structure  of  the  field  as  it  disregards  the  overlapping  and  interaction  of 
certain  wave  types.  However,  with  respect  to  the  energy  input,  we  can 
easily  tolerate  these  adulterations.  To  demonstrate  the  influence  of  the 
sediment  layer  we  keep  the  source  position  fixed  in  the  middle  of  a  60  m 
water  column  and  decrease  the  layer  thickness  stepwise.  In  Fig.  14  we  have 
plotted  the  transmission  loss  against  range  in  the  acoustic  (dashed  curves) 
and  the  seismic  (solid  curves)  fields  for  the  cases  of  two  enormous  sand 
layers  of  120  m  and  60  m  thickness.  In  both  cases  the  acoustic  propagation 
is  quite  bad  and  the  seismic  one  even  worse,  but  we  realize  at  once  that  a 
decrease  in  the  layer  thickness  slightly  improves  the  conditions  for 
seismic  sensing  and  worsens  those  for  acoustic  detection.  A  more  detailed 
analysis  reveals  that  the  acoustic  f^eld  comprises  merely  one  mode-like 
structure,  while  there  are  different  weak  contributions  to  the  seismic 
field.  Among  the  latter,  the  Scholte  wave  at  the  water/sand  interface 
(csch~600  m/s)  can  be  identified  clearly,  but  due  to  the  high  absorption 
in  the  sediment  it  is  of  no  importance. 


RANGE (km) 

FIG.  14  RANGE  DEPENDENCE  OF  THE  ACOUSTIC  (DASHED)  AND  SEISMIC  (SOLID) 
TRANSMISSION  LOSS  FOR  A  10  Hz  SOURCE  OVER  DIFFERENT  THICK 
SAND  LAYERS 
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Omitting  some  intermediate  steps,  which  confirm  the  above-mentioned  trend, 
we  pass  immediately  to  three  relatively-thin  sand  layers  of  10  m,  5  m,  and 
0.5  m  thicknesses.  The  pertinent  curves  in  Fig.  15  demonstrate  impressively 
that  seismic  sensing  can  offer  not  only  an  alternative  but,  on  the 
contrary,  may  repres<?r\  a  superior  detection  method.  The  remarkable  ranges 
achieved  by  seismic  propagation  are  solely  a  consequence  of  a  modified  or 
generalized  Scholte  wave  that  is  guided  along  the  "sand-coated"  water/rock 
interface  (csc^  «  1150  m/s). 

It  should  be  mentioned  that  there  is  no  direct  contribution  to  the  seismic 
field  from  the  side  of  the  sand/rock  interface,  because  our  choice  of 
parameters  excludes  the  existence  of  a  "pure"  Stoneley  wave  at  this 
boundary  (see  Sect.  1.1). 


RANGE (km) 

FIG.  15  RANGE  DEPENDENCE  OF  THE  ACOUSTIC  (DASHED)  AND  SEISMIC  (SOLID) 
TRANSMISSION  LOSS  FOR  A  10  HZ  SOURCE  OVER  DIFFERENT  THIN 
SAND  LAYERS 
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The  excitation  of  the  modified  Scholte  wave,  and  thus  the  good  seismic 
propagation,  would  be  improved  if  the  source  position  were  shifted  towards 
the  ocean  floor.  For  that  reason  Fig.  16  compares  the  seismic  transmission 
loss  of  the  preceding  plot  with  the  corresponding  curves  calculated  for  a 
bottom-mounted  transmitter  SQ  of  the  same  source  level.  In  all  three  cases 
a  decrease  of  the  loss  by  about  10  dB  confirms  the  expected  better  coupling 
of  source  and  sea  bed. 


0  10  20 

RANGE (km) 

FIG.  16  RANGE  DEPENDENCE  OF  THE  SEISMIC  TRANSMISSION  LOSS  FOR  A  10  Hz 
SOURCE  OVER  AND  ON  A  THIN  SAND  LAVER 


Finally,  to  show  that  the  selected  water  depth  H-j  of  merely  60  m  does  not 
imply  a  serious  limitation  to  our  conclusions,  we  returned  to  the  inter¬ 
mediate  example  of  Fig.  15  (source,  SQ,  30  m  above  a  5  m  sand  layer)  and 
calculated  the  transmission  loss  for  water  columns  of  80  m,  100  m,  and  120  m. 
From  the  resulting  curves  of  Fig.  17  we  see  that  the  poor  acoustic 
propagation  improves  only  slightly  with  increasing  water  depth,  while  the 
relatively  good  seismic  propagation  is  not  affected. 
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To  give  an  estimation  of  the  optimum  frequency  range  for  seismic  detection 
on  top  of  such  a  "sand-coated"  rock  bottom  (intermediate  case  of  Fig.  17) 
the  frequency  was  varied  instead  of  the  range  (receiver  distance  10  km). 


RANGE (km) 

FIG.  17  RANGE  DEPENDENCE  OF  THE  ACOUSTIC  (DASHED)  AND  SEISMIC  (SOLID) 

TRANSMISSION  LOSS  FOR  A  10  HZ  SOURCE  IN  DIFFERENT  WATER  COLUMNS 
OVER  A  THIN  SAND  LAYER 

Figure  18,  which  stems  from  another  study  [100),  shows  the  pertinent 
spectrum  levels  of  the  acoustic  and  the  seismic  transmission  losses,  if 
the  frequency  is  decreased  below  30  Hz  the  acoustic  propagation  loss 
increases  drastically  by  about  40  dB,  with  a  relatively  broad  maximum 
between  20  and  10  Hz.  At  about  5  kHz  the  acoustic  propagation  again 
reaches  its  former  "high-frequency"  quality  in  the  form  of  a  relatively 
sharp  loss-minimum,  but  then  it  is  more  or  less  completely  extinguished 
below  3  Hz.  This  narrow-band  improvement  around  5  Hz  is  due  to  the 
incorporation  of  the  sand  layer  into  the  acoustic  duct. 
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FIG.  18  FREQUENCY  DEPENDENCE  OF  THE  ACOUSTIC  (DASHED)  AND  SEISMIC  (SOLID) 
TRANSMISSION  LOSS  AT  A  FIXED  RANGE  OF  10  km  FOR  A  10  HZ  SOURCE 
OVER  A  THIN  SAND  LAYER 


On  the  other  hand,  the  initially  very  high  seismic  propagation  loss  falls 
below  its  acoustic  counterpart  at  about  20  Hz,  then  drops  for  another 
50  dB,  and  finally  passes  a  relative  flat  minimum  between  7  and  2  Hz.  Thus 
for  each  tonal  below  20  Hz  the  seismic  propagation  is  much  better  than  its 
acoustic  counterpart. 


An  additional  support  for  the  better  seismic  sensing  (better  signal-to- 
noise  ratio)  may  result  from  a  typical  minimum  of  the  ambient  noise 
spectrum  levels  in  the  10  to  20  Hz  band  (see  Sect.  2.3). 


SACLANTCEN  SR-42 


CONCLUSIONS  AND  PROPOSALS  FOR  FUTURE  RESEARCH  PROGRAMMES 

A  substantial  part  of  the  world's  oceans  can  be  classified  as  shallow  water 
(e.g.  water  depth  of  200  m  or  less),  and,  what  is  more  important,  many  of 
these  areas  (the  continental  shelfs  and  most  of  the  straits)  are  strategi¬ 
cally  important.  It  is  in  just  these  areas  that  the  approved  and  highly- 
developed  sonar  techniques  for  long-  and  medium-range  detection  lose  their 
accuracy  or  even  fail. 

When  propagating  in  shallow  water,  or  from  deep  to  shallow  water  over  a 
sloping  bottom,  the  acoustic  energy  of  "water-waves"  becomes  rapidly 
attenuated  through  repeated  interactions  with  surface  and  bottom,  similarly 
the  coherence  of  the  pertinent  sound  field  is  also  lost.  On  top  of  that, 
all  very  shallow  regions,  and  thus  areas  close  to  shore,  are  impenetrable 
for  infrasonic  energy  (approximately  the  1  to  30  Hz  frequency  band)  due  to 
the  existence  of  an  absolute  low-frequency  cut-off  for  the  water  wave 
(meaning  the  one  for  the  first  or  lowest  propagation  mode  of  this  field). 

Thus  the  conventional  acoustic  surveillance  techniques  suffer  from 
shortcomings  (on  the  part  of  both  hardware  and  software)  and  ask  for 
alternatives.  Besides  different  non-acoustic  methods  (based  on  optical, 
magnetic,  or  chemical  phenomena)  a  closely-related  alternative  may  be 
offered  by  that  fraction  of  the  acoustic  energy  that  becomes  coupled  with 
the  bottom  and  is  thus  transmitted  in  the  form  of  seismic  waves  along  and 
in  the  sea  floor.  Although  seismology  has  been  an  extensive  field  of 
research  for  many  years  this  coupling  mechanism  is  still  poorly  understood 
(especially  with  respect  to  interface  waves),  not  to  mention  the  very 
confusing  variety  and  interaction  of  different  propagation  modes.  That  is 
why  we  tried  to  give  a  comprehensive  qualitative  description  of  the 
underlaying  physics  in  Ch.  1.  Despite  the  use  of  directional  sensors  and 
highly-sophisticated  signal-processing  techniques  in  seismology,  the 
related  and  very  crucial  problems  of  signal  strength  and  coherence  have 
been  largely  neglected  due  to  the  fact  that  the  questions  of  energetics  are 
of  minor  interest  in  prospecting  and  earthquake  monitoring.  Nevertheless 
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we  are  now  able  to  extract  numerical  results  under  strongly  simplified 
conditions  with  the  help  of  modern  computer  programs,  as  has  been  shown  in 
Ch.  3. 

Unfortunately  these  latter  possibilities  are  not  only  restricted  by  the 
natural  limitations  of  the  existing  models  but  also,  and  even  more,  by  our 
lack  of  knowledge  concerning  the  environmental  input  data. 

The  large  amount  of  published  data  on  the  lithosphere  below  the  oceans  can 
merely  indicate  the  limits  and  trends  for  modelling  work  and  predictions. 
That  is  to  say,  the  seismic  parameters  involved  are  not  only  more  numerous 
but  are  also  affected  by  much  stronger  variations  with  depth  and  range  due 
to  (non-regular)  layering  and  inhomogeneities.  Especially  with  respect  to 
reliable  data  on  shear  phenomena  in  sediments,  we  are  still  largely  groping 
in  the  dark.  Concerning  the  seismic  ambient-noise  levels  in  shallow  water, 
the  starting  conditions  are  even  worse.  The  very  few  existing  measurements 
with  vertical  geophones  differ  considerably  and  do  not  permit  us  to  draw 
any  conclusions  with  regard  to  the  more  interesting  horizontal  unrest  of 
the  sea  floor  (including  the  questions  of  directionality).  Thus  the  only 
information  on  favourable  signal-to-noise  ratios  of  bottom-mounted 
horizontal  geophones  (exceeding  those  of  nearby  mounted  hydrophones  by 
about  5  to  20  dB)  stems  from  a  few  non-coordinated  field  trials  in  the  past 
(see  e.g.  [182-184]). 

Likewise,  the  coherence  characteristics  of  the  wave  types  in  the  marginal 
portion  of  the  seismic  field  at  or  just  below  the  sea  floor  (mainly 
interface  and  trapped  shear  waves)  are  widely  unknown  but  it  has  to  be 
expected  that  they  are  strongly  deteriorated  by  lateral  inhomogeneities 
[198]. 

In  Ch.  2  we  tried  to  shed  some  light  on  the  historical  and  technical 
reasons  for  all  these  short-comings  and  contradictory  results  in  former 
studies.  As  a  result,  we  can  define  the  outline  of  future  research 
activities  by  listing  the  necessary  experimental  and  theoretical  tasks  in 
order  of  importance. 
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More  careful  and  systematic  measurements  should  be  performed  with  OBS- 
stations  in  shallow  water  to  identify  and  explore: 

a.  The  excitation  and  propagation  mechanisms  of  the  different 
interface  phenomena  (pertinent  levels,  phase  and  group  velocities,  losses, 
and  directionality). 

b.  The  levels  and  directional  characteristics  of  acoustic  and  seismic 
ambient  noise  in  the  infrasonic  frequency  range  (about  1  to  30  Hz). 

c.  The  achievable  improvements  of  the  signal-to-noise  ratio  by  using 
horizontal  geophones  instead  of  or  together  with  omnidirectional  hydrophones. 

d.  The  bottom  characteristics  of  the  top  layers  in  general  and  the 
shear  damping  in  particular. 

e.  The  deterioration  of  signal  coherence  by  lateral  inhomogeneities. 

These  experimental  studies  should  be  accompanied  by  a  lot  of  interacting 
modelling  work  to  facilitate  the  interpretation  of  field  data  and  to 
support  the  selection  of  suitabk  test  areas. 
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APPENDIX  A 

POST-1978  WORK  AT  SACLANTCEN 


INTRODUCTION 


The  main  text  was  drafted  in  late  1978  but  other  priorities  delayed  its 
publication  until  the  end  of  1980.  Thus  its  references  are  to  the  open 
literature  published  before  1979  and  the  comments  made  in  Sect.  2.4  about 
SACLANTCEN's  unsuccessful  field  work  are  valid  only  for  the  early  stages  of 
SACLANTCEN' s  project  on  bottom  interface  waves. 

In  late  1978  SACLANTCEN's  tri-axial  sensor-package  and  its  radio-buoy  were 
completely  re-designed  and  the  first  version  of  a  new  system  was  built  in 
early  1979.  The  following  pages  briefly  describe  these  more  advanced  tools 
and  present  a  few  typical  results  from  our  1979  sea  trials. 


A.  1  New  Equipment 

Besides  a  change  to  an  all-digital  technique  the  main  break-through  was 
obtained  by  using  a  new  type  of  active  and  lightweight  1  Hz  geophones 
(Teledyne  Geotech  S-500).  These  sensors  do  not  need  blocking  for  trans¬ 
portation  and,  because  they  can  be  applied  in  any  orientation,  do  not  have 
to  be  precisely  levelled  for  proper  operation.  The  seismic  sensor  set  is 
completed  by  a  variable-depth  hydrophone  mounted  outside  the  container  of 
the  ocean-bottom  seismometer  (OBS)  or  floating  above  it.  Figure  A1  shows 
how  this  sensor  package  is  installed  in  shallow  water  together  with  its 
radio  buoy. 


FIG.  A1  INSTALLATION  OF  THE  SENSOR  PACKAGE  ON  THE  SEA  FLOOR 
AND  MOORING  OF  ITS  RADIO  BUOY  IN  SHALLOW  WATER 
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Figure  A2(a)  presents  a  simplified  block-diagram  of  the  pertinent 
electronics,  which  are  similar  to  those  described  by  Blackinton  and  Odegard 
<A.1>. 

The  output  of  each  of  the  four  basic  sensors  is  preamplified,  low-pass 
filtered  (100  Hz),  and  fed  direct  and  via  three  additional  amplifier-stages 
of  18  dB  each  into  a  4  x  4  channel  multiplexer.  Each  of  the  four  signal 
levels  is  scanned  with  a  sampling  frequency  of  600  Hz,  and  the  optimally- 
amplified  signal  is  transmitted  to  a  fast  12-bit  A/D  converter.  The 
resulting  12-bit  mantissa  is  then  combined  with  the  pertinent  2-bit 
exponent,  with  a  parity  bit,  and  with  another  independent  bit  to  a  basic 
data  word  of  16  bits.  After  three  scanning  cycles  of  the  basic  sensors  we 
thus  have  available  twelve  additional  bits  for  the  transmission  of 
important  parameters  and  control  data.  Eleven  of  those  are  used  to  form  a 
secondary  data  word,  which  is  cyclically  assigned  to  one  of  the  four 
auxiliary  sensors  (compass,  tilt  X,  tilt  V  and  temperature)  and  one  of  the 
four  electronic  checkpoints..  Thus,  every  24  basic  scanning  cycles  we  can 
not  only  monitor  the  actual  position  of  the  OBS,  the  nearby  water 
temperature,  and  the  power  supply  of  four  important  networks,  but  also  have 
at  our  disposal  an  additional  set  of  eight  bits  to  use  for  the  OBS 
station's  identification  code.  These  digital  data-sequences  are 
transferred  to  the  surface  buoy  via  a  coaxial  cable  and  from  there  to  the 
receiving  ship  by  a  FM-modulated  radio-transmitter  (170  MHz,  max  15  W)  in 
the  surface  buoy.  This  buoy  also  houses  the  rechargeable  batteries  that 
power  the  complete  system  and  a  radio  receiver  for  the  ON-OFF  commands  sent 
from  the  receiving  platform.  Every  time  the  OBS  is  switched  on,  a 
calibration  signal  is  automatically  applied  to  each  geophone  calibration 
coil  for  about  one  to  two  minutes.  One  of  the  most  prominent  features  of 
this  digital  OBS  is  given  by  the  very  high  dynamic  range  of  120  dB  (66  dB 
from  the  12-bit  mantissa  and  54  dB  from  the  2-bit  exponent),  which  enables 
us  to  cope  simultaneously  with  a  low-level  seismic  background  and 
high-level  deterministic  signals  from  nearby  CW-sources  or  explosions. 

Figure  A2(b)  shows  the  corresponding  block-diagram  of  the  data-acquisition 
and  pre-processing  facilities  on  board.  Here  we  are  using  a  Hewlett- 
Packard  21  MX  computer  with  disc  unit,  which  gives  us  the  possibility  of 
getting  a  printout  or  a  plot  of  all  essential  data  channels  a  few  seconds 
after  recording  an  event. 
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FIG.  A2  ELECTRONIC  BLOCK-DIAGRAM  FOR: 

(a)  the  OBS  and  its  radio-buoy, 

(b)  the  onboard  receiving,  recording,  controlling 

and  pre-processing  systems 
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A. 2  Experimental  Results 

Figure  A3  indicates  two  positions  off  the  Italian  coast  where  the 
instruments  were  deployed  during  a  cruise  in  1979.  The  sensor  station  was 
always  oriented  in  such  a  way  that  the  directions  of  the  planned  acoustic 
runs  (at  constant  water  depth  or  a  sloping  sea  floor)  coincided  with  the 
axis  of  one  of  the  two  horizontal  geophones.  As  sound  sources  we  used 
small  TNT  charges,  which  were  usually  placed  on  the  sea  floor  and  fired 
electrical ly. 


FIG.  A3  THE  OPERATIONAL  AREA  OFF  THE  VERSILIAN  COAST  (ITALY) 
WITH  THE  TWO  CHOSEN  OBS-POSITIONS 
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To  discuss  some  of  the  typical  features  of  the  interface  waves,  we  will 
look  first  at  an  event  from  Run  3  along  the  20  m  depth  profile  off 
Viareggio,  where  the  hydrophone  was  floating  in  the  middle  of  the  water 
column.  Figure  A4  presents  a  lineprinter  plot  of  the  four  time-series 
generated  by  a  charge  of  180  g  TNT  at  a  distance  of  about  1.3  km.  The 
lower  three  traces  display  the  measured  particle  velocities,  u,  v  and  w, 
while  the  upper  trace  shows  the  pressure  history  recorded  by  the 
hydrophone.  The  latter  channel  was  always  subject  to  an  additional 
pre-amplification  of  30  dB.  Due  to  this  scaling  factor  we  usually  had  to 
clip  the  first  arrival,  which  was  formed  by  the  unavoidable,  high-frequency 
water  wave.  On  this  and  all  following  plots  the  time-series  were  passed 
through  a  10  Hz  low-pass  filter  to  stay  below  the  cut-off  frequency  of  the 
water  duct. 

The  interface  wave  we  are  looking  for  arrives  somewhat  later  than  the  water 
wave.  As  expected  from  our  geometrical  consideration,  this  interface  wave 
is  detected  only  by  the  vertical  geophone  GeoZ  and  the  radial  one  GeoX, 
which  is  horizontal  and  oriented  parallel  to  the  direction  of  propagation. 
The  ripples  on  the  output  of  the  transverse  sensor  GeoY  are  probably  due  to 
imperfections  of  the  wave  guide  and  lateral  inhomogeneities. 


FIG.  A4  ^INE-PRINTER  PLOT  OF  THE  SIGNALS  FROM  THE  FOUR  BASIC 
SENSORS  (180  %  TNT  FIRED  AT  1.3  km  DISTANCE) 
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The  wavelet  itself  consists  of  a  long-lasting,  narrow-band  signal  that 
shows  the  expected  normal  dispersion,  indicating  that  the  higher 
frequencies  are  left  behind.  Especially  at  the  output  of  the  radial  sensor 
GeoX  we  realize  that  the  beginning  of  the  main  signal  is  concealed  by 
another  superimposed  wavelet  that  has  relatively-high  amplitudes  and  an 
extremely  small  bandwidth.  In  agreement  with  the  theoretical  predictions, 
the  hydrophone  also  senses  the  low-frequency  pressure  variations  created  by 
the  interface  wave. 

As  a  final  check  for  identification  we  have  displayed  the  corresponding 
velocity  vector  in  the  radial/vertical  X-Z  plane  using  the  radial  velocity 
u  as  abscissa  and  the  vertical  velocity  w  as  ordinate.  At  the  top  of  Fig. 
A5  we  have  therefore  again  plotted  the  appropriate  sections  (windows)  of 
both  signals  and  below  we  have  displayed  the  resulting  hodographs  for  14 
successive  one-second  windows. 

De''  to  the  above-mentioned  superimposed  wavelet,  the  first  four  windows  A-D 
are  still  quite  irregular  and  curve  B  especially  seems  to  be  complicated  by 
the  superimposition  of  a  horizontal  ellipse.  This  latter  feature  could 
indicate  the  coexistence  of  a  higher  mode.  Nevertheless,  the  following  ten 
hodographs  (E  to  N)  form  very  regular  vertical  ellipses  that  are 
circumscribed  clockwise  or  prograde.  Towards  the  tail  of  the  signals  these 
ellipses  fade  away  and  merge  into  the  background  noise. 
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FIG.  A5  RADIAL  AND  VERTICAL  PARTICLE  VELOCITY  OF  THE  INTERFACE 
WAVELET  IN  FIG.  A4  WITH  THE  RESULTING  HODOGRAPHS 
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After  this  unambiguous  identification  of  the  wave  type  we  would  like  to 
discuss  more  carefully  some  of  its  propagation  characteristics.  For  this 
purpose  we  have  stacked  in  Fig.  A6(a)  the  output  of  the  vertical  geophone 
GeoZ  for  Run  2,  using  the  same  scale  for  all  signals  but  increasing  the 
actual  charge  size  with  growing  distance  according  to  the  indication  on  the 
right-hand  side. 
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FIG.  A6  (a)  STACKING  OF  THE  VERTICAL  PARTICLE-VELOCITY  FOR 
AN  ACOUSTIC  RUN  OVER  A  LAYERED  SEDIMENT-BOTTOM 


As  reference  quantity,  the  first  dotted  line  marks  the  onset  of  the 
unavoidable  water  wave  propagating  with  a  velocity  of  1536  m/s.  The  first 
continuous  line  indicates  the  front  of  the  seismic  signal,  which  is  roughly 
characterized  by  1.5  to  2  Hz  oscillations  and  a  velocity  of  about  240  m/s, 
while  the  second  continuous  line  marks  the  tail  of  the  signal  with  4 
to  5  Hz  oscillations  and  a  velocity  of  75  m/s.  In  our  plot  this  latter 
tail  of  the  seismic  wavelet  disappears  at  distances  of  about  1  km  due  to 
the  frequency-dependence  of  the  attenuation;  this  dependence  may  be  linear, 
as  is  known  from  P-  and  S-wave  propagation  at  higher  frequencies.  The 
attenuation  coefficients  for  this  and  several  other  data  sets  will  De 
calculated  properly  and  reported  separately. 

The  seismic  wave-guide  under  consideration  thus  offers  optimal  propagation 
conditions  for  frequencies  below  3  Hz,  The  second  dotted  line  emphasizes 
the  above-mentioned  superimposed  wavelet,  which  has  a  propagation  velocity 
of  185  m/s. 
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Figure  A6(b)  shows  the  stacking  of  the  corresponding  output  of  the  radial 
geophone  GeoY.  Here  the  onset  of  the  interface  wave  is  more  hidden  by  some 
refracted  arrivals  but  the  superimposed  wavelet  is  much  more  pronounced  and 
reveals  its  very  narrow  frequency  band  of  2.5  to  3  Hz. 

Modelling  this  special  wave-propagation  phenomenon  with  the  Centre's 
FFP-Program  showed  that  the  theory  indeed  predicts  seismograms  composed  of 
a  long-lasting  dispersed  wave  train  (zeroth  mode)  and  a  super-imposed, 
pulse-like,  more-or-less  monochromatic  wavelet  (first  mode).  A  more 
detailed  discussion  of  our  measurements  and  their  theoretical  aspects  was 
presented  <A.2>  at  the  conference  on  Bottom-Interacting  Ocean  Acoustics 
Conference  held  at  SACLANTCEN  in  June  1980.  In  a  second  paper  <A.3>  at  the 
same  conference  we  demonstrated  that  in  shallow-water  areas  ambient  and 
ship- induced  infrasonic  noise  is  transmitted  in  the  form  of  Scholte  waves 
and  thus  offers  the  possibility  of  determining  the  bearing  of  the  source. 


FIG.  A6  (b)  ST  'KING  OF  THE  RADIAL  PARTICLE-VELOCITY  FOR  THE 
SAME  ACOUSTIC  RUN  AS  IN  A6(a) 


SACLANTCEN  SR-42 


i 


i 


e 


iz.  »^S 


i 

\l 

r 

.1 

j 

t 

\ 


REFERENCES 

A.  1  BLACKINGTON,  J.G.  and  ODEGARD,  M.E.  An  ocean-bottom  seismograph  using 
digital  telemetry  and  floating-point  conversion.  IEEE  Trans. 
Geoscience  Electronics  GE-15:  74-82  (1977). 

A. 2  RAUCH,  D.  Experimental  and  theoretical  studies  of  seismic  interface 
waves  in  coastal  waters.  In:  KUPERMAN,  W.A.  and  JENSEN,  F.B.,  eds. 
Bottom- interacting  ocean  acoustics  conference,  June  9-13,  1980. 

SACLANT  ASW  Research  Centre,  La  Spezia,  Italy.  (To  be  published  by 
Plenum  Press). 

A. 3  SCHMALFELDT,  B.  and  RAUCH,  D.  Ambient  and  ship-induced  low-frequency 
noise  in  shallow  water.  In:  KUPERMAN,  W.A.  and  JENSEN,  F.B.,  eds. 
Bottom- interacting  ocean  acoustics  conference,  June  9-13,  1980. 

SACLANT  ASW  Research  Centre,  La  Spezia,  Italy.  (To  be  published  by 
Plenum  Press). 


103 


f 


jv 


DATE 

FILMED 

/ 


